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A general structural theory of fields which arise from variational statements 
in a non-metric but affinely connected space without torsion is shown to be a 
natural setting for an intrinsic geometric description of electromagnetic, gravi- 
tational, and matter fields. No specific intrinsic geometry beyond that of an 
affinely connected space is assumed in the derivation of the field equations, 
and no particular transformation properties are assumed for the Lagrangean 
function. The field equations which result are thus not necessarily tensorial in 
nature, although equations of Maxwellian form follow in a natural way together 
with definite constitutive relations for the charge-current density and the charge- 
current potential. An electromagnetic field tensor is obtained even if the 
field equations be nontensorial and is shown to be twice the skewsymmetric 
part of the third contraction of the affine curvature tensor. This field tensor 
is fundamental in the development, and it is shown that its vanishing is a necessary 
condition for the affinely connected space to reduce to a metric space. Under 
an explicit definition of ‘“‘emptiness” the general affinely connected space is 
shown to reduce, for a particular Lagrangean function, to a metric space in 
which EINSTEIN’s equations for the vacuum gravitational field are satisfied. 
For a general class of Lagrangean functions the space which is defined as “‘semi- 
empty” is shown to be a manifold with a Weyl geometry; and for a particular 
Lagrangean function of this class the field equations describe classic vacuum 
electrodynamics. It is then shown that the classic vacuum electrodynamic and 
gravitational fields are a simple separation of effects. 
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1. Notation 

A point set D contained in a space & will be referred to as a domain if it is 
an open, connected, non-empty point set. We denote by D* the closure of D. 
By the boundary of D* shall be meant the set D* ©D, where © stands for the 
set-theoretic difference. 

Partial differentiation will be denoted by a comma, and covariant differen- 
tiation relative to the affinity L%, will be denoted by a semi-colon. When a 
partial derivative is to be taken with respect to an element of a collection of 
functions defined over D or D*, it will be denoted by a comma followed by the 
symbol for the particular element as a subscript (H,,=06H/@¥). 

We adopt the following symbolic representations for volume and surface 
differentials: 


N,dS2dddw?..dv 4d? .... 


The symmetric and skewsymmetric parts of a tensor T,, will be denoted 
by Tis) and T,,¢) respectively. 


2. The Affine Curvature Tensor and Its Contractions 


Let & be a four-dimensional Hausdorff space into which a structure is intro- 
duced by requiring that at every point of & there be defined a symmetric affine 
connection L%,. The affine curvature tensor K’,,,,, associated with L%,,, is 


defined by 


Ly a tania as Laem © ee oo ee (2.1) 
from which 
AN su Mus = BO  A® (2.2) 


where A* are the components of a contravariant tensor of rank one and weight 


zero. We may build three new tensors from the tensor K”,,,, by contraction. 


aX fl 
The first contraction: 


def ;> 
Ag = IO eo pe ai (2.3) 
The second contraction : 
def r x ” Y 7 
Bei a9 Oa = Len i ci Den as reas L'bx rT Le Lu (2.4) 
The third contraction: 
def y> Y ; 
Ce ees Ioan re De F Doe = Lig die — TEP) Day: (2.5) 


. The three tensors defined by equations (2.3) through (2.5) are not algebraically 
independent. Specifically, since the affine connection of & is assumed to be 
symmetric, the affine curvature tensor admits the following complete set of 
algebraic identities!: 


a | 
Kun) = 9, ES tees 


1 THOMAs, T. Y.: The Differential Invariants of Generalized Spaces, p. 132. Cam- 
bridge: Cambridge University Press 1934. 
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Substituting equations (2.3) through (2.5) into these identities, we have 
Beenie BAG OCs a. 


Thus, C,,, is the only independent contraction of the affine curvature tensor, 
the other two contractions being obtained from C, . by the above equations. 
Rather than work with the tensor C,,, directly, we decompose C,,,, into its skew- 
symmetric and symmetric parts, as defined by 


Agu= eP BG fe sf dens py — Le, (2.6) 
and 


def € € é € € 
Be = Cap) = 3 (L caer ena) —_ als tg Ba Leap (2.7) 


respectively, and work with these tensors separately. As will be seen in Section 4, 
this decomposition of C,,,, admits a simple and direct separation of electromagnetic 
and non-electromagnetic effects. 

The vanishing or non-vanishing of the tensor A,, determines certain in- 
trinsic properties of the space & which will be required in succeeding sections. 


Theorem 1. A covariant constant scalar density? field 1s possible in & if and 
only if A,g=0. 

Proof. By definition, a density C(#) forms a covariant constant field in & 
if and only if 

Caste TOS 
at all points of &. From the above equation we obtain 
(nC) ,=L*.,. 

These equations are integrable for InC, and hence for C, if and only if L*;,, ,;=0, 
which implies A,,,,=0 by equation (2.6). 

Theorem 2. A covariant constant quadrivector field is possible in & 1f and 
only if Ayg=0. 

Proof. The proof, using the integrability theorem for mixed systems of partial 


differential equations, is like that of Theorem 4. 


An affinely connected space is said to be a metric space if there is defined 
at all points a tensor g, satisfying the following postulates: 


(i) Stag) = 9, 
(il) g = det (g.) +0, 
(iii) Lagiy = 04, (2.10) 


in which case the fundamental quadratic form is given by gy, dx” dx? 


2 A quantity C(a) is said to form a scalar density if, under a transformation of 
x of class C1, it transforms according to the law ’C = [det (‘a ,,)]1C. This definition 
is included since the term “‘scalar density’ is used by different authors to mean 
several different things. 

8 SCHOUTEN, J. A.: Ricci-Calculus, pp. 154—155. Berlin-Gottingen-Heidelberg: 
Springer 1954. 
4 SCHOUTEN: op. cit., p. 132. 
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Theorem 3. The affinely connected space & cannot reduce to a metric space 
tf Ags. 

Proof. It can be shown that a necessary condition for the existence of a 
tensor field g,, satisfying the postulates of a metric space is expressed by the 


equations ® : 
4 Bop Kaya t SueK pys=9- 


Multiplying by g’’ and contracting with respect to (A,«) yields 2K%,,5=0. This 
implies, by equation (2.3), that A,;=0, from which the theorem follows. 
Let there be defined over an affinely connected space a tensor field g,, satis- 


fying the following postulates: 


. Poy ee (2.114) 
(ii) g = det (81) +0, Cou 
(iii) Bapsy = — Oyap igs 

where (i) and (iii) imply ieee On (2.14) 


Under these postulates it can be shown® that a symmetric affine connection L*,; 
is uniquely determined and has the form 


Pee Ghee” (Giger Cre ean) (2.15) 
where /™,, are the Christoffel symbols based on g,¢, 7.e., 
"k= hier Sas Cane (2.16) 
and a 
Ia = (In |g1),2- (2.17) 


From equation (2.15) we have 
(wa) = 38° (Qfuate + Ole! al) 
upon using equation (2.14). Thus L%,, is symmetric if and only if 
Oluale + rately = 0- 
Rearranging the terms in the last equation gives 
Qutae] = Qatnel 


which is identically satisfied since both terms vanish by equation (2.14). Thus 
equation (2.15) defines a symmetric affine connection. Set 


Le aes 5 (Orne Qron Qo a) , (2.18) 

and denote covariant differentiation relative to the connection I, by V >: thent 
—_ ¢ C 

Rua=JIuat Veg? IB 1 ee Teag — $V 8% Teno + 

get 9° * Tine Tameree ona Trae, 


Mee 


(2.19) 


a S 
AT, ease 2(g® Liu \tol), 21 (2.20) 


° THomas, T. Y.: of. cit., pp. 218—219. (The B® fol 1 

en ae ( ayo Of THomas is the same as the 
§ SCHOUTEN, J. A.: op. cit., p. 132, equation 3.5. 
” SCHOUTEN, J. A.: op. cit., p. 141. 
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where J,, is the third contraction of the affine curvature tensor associated with 
the connection 1%... 


Let OL’, be an arbitrary variation of L%, which is constrained by the con- 
dition that the symmetry of L%,,, shall be preserved under the process of variation, 
v.e., OL%,—=dL*,,. We may easily show, by direct computation, that the follow- 


It holds: 
ing result holds OK a (OL) 410s), 


we? 
from which we obtain 


0Axg = (Lp) ,0— (OL%0);8, (2.21) 
OL Olen AOL 8) a (Olewa). o: (2.22) 


3. Variational Considerations and the Affine Field Equations 

Let Qi.) be a collection of functions, which we shall call the exterior matter 
fields; their transformation properties are left unspecified. Let Y(A,,, Ryg, 
L534, Qe), Qo),e. *) be a function whose transformation properties are unspecified 
and which is of class C?in its arguments; this function will be called the Lagrangean. 
Although we could also include L%,,, as an explicit argument of ¥ in addition 
to its implicit occurrence through A,, and R,,, this additional complexity is 
found unnecessary. We wish to determine L%,, and Q,,) such that 


6f GdV=0 (3.1) 
D* 


for all variations 6L%,, and 0Q,,) which vanish on the boundary of D* and are 
such that dLfg.,=0. 


Define the quantities /*° and H*? by 
LP = Lo, 6 (3.2) 
HP =f 4, (3.3) 


which are symmetric and skewsymmetric respectively, due to the respective 
symmetry and skewsymmetry of R,, and A,,. Using equations (2.21) and (2.22), 
we obtain, by direct computation, 


es [Zoe (LF oO, ¢) o ce] OQ dV + 


$f [HOM OL ea Ne — 9M) + (oy,0) 900M + 
+ 1°82 61% 2N, +4 06L%.N, — 6146 N,)|dS + 
ety J ee Oy i Hats = EPI, pee abe 


D* 


PPD — 1° On 2H”* 0, + O19, 2 | OL, dV 


where Or, f is evaluated holding A,, and R,, constant. This result can be 
simplified significantly by interpreting I*® as if it were a contravariant tensor 
of rank two and weight one and then applying the operation of covariant differ- 
entiation with respect to L%,, 1.¢., 


ited cee Lg Det Lie lt — bi, Leh. 
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It must be noted that [%° is not necessarily a tensor density; it will be one only 
if is a scalar density, but this is not assumed at present. With this operational 


agreement, we have 
bf GIV=J [Lo — (Lo,e),t] 64 V + 
D* D* 


+ f [01% 2+ Pg — BT. + 26,H" | OL, dV + 
D* 


; (3.4) 
+ f [H*?(6L%—N, — OL %sq Np) + 


D*¥OD 


ae T° (£6 Leg Ng =i z Oleea Ne y OL, g N,) ate F 0,t 6 Qe) N;| ds. 


Requiring the variations to vanish on the boundary of D* and 6L%%,,;=0, we 
obtain the following necessary conditions for (3.1) to hold: 


DP nes oy TP* , — Gl ae 2078, aL 208 H’*,+ 20H" = 0, (3.5) 
Z 0) (LF o,¢),¢ = 9- (3.6) 


Equations (3.5) and (3.6) are 40+ (g) equations for the determination of the 
40+ (0) field variables’ L%, and Q,,), and are the basic affine field equations 
of this theory. 


Theorem 4. Equation (3.5) is equivalent to the equations 
H*.= 10 BI ie 244, 2), (3. 
(PY dp — BLP Oe — BI Or) = (On 18, Bon Ong — 1%) Se 8) 


Proof. Contracting equation (3.5) with respect either to («, ) or to (a, y) 
gives equation (3.7). Eliminating H%® between equations (3.5) and (3.7) yields 
equation (3.8). 

Note. If we had not assumed that the torsion tensor of & vanishes, then 
contracting the equation which would replace equation (3.5) with respect to 
(x, 8B) and («,y) would give two different right-hand sides to equation (3.7). 
Thus the resulting equation (3.5) would be replaced by the system composed 
of equations (3.7), (3.8), and the equation which would result from equating 
the two resulting right-hand sides of equation (3.7). Such conditions would 
yield an additional four equations. In this case we should then have 68+ (@) 
equations to determine 64+ (0) variables. 


4. Equations of Maxwellian Form as Natural Consequences of the 
Affine Field Equations 


In tensorial form, MAxXWELL’s equations for a linear medium are 
Kapy=9, fap, =9: 
Vo a0 Le =; 

op é 

Fi er 


22? and ee are known functions of L%,-, Qo) and their derivatives once the 
Lagrangean function ¥ has been given (cf. equations (3.2) and (3.3)). In this sense 
the Lagrangean function plays the role of a constitutive function. 
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where G%8”® is a tensor representing the properties of the medium. More generally, 
we shall say that a system is of Maxwellian form if it may be written as 


hap = 0, hap.) = 9; 
REM S10 poleRreye st, (4.1) 


Tel ag A ges +) typ? 


We shall refer to J(/,,,...) as the constitutive function for the Maxwellian field. 
The tensor A,, satisfies the conditions 


Aggy = 0, Atas,yj]= 9, 


as is evident from the defining equation (3.6). We have seen that upon con- 
tracting the affine field equations we obtain 


H®* = io (31%*. pit. 261°, £) 
(cf. Theorem 4). Setting 
Sg Bl te 202, 2), (4.2) 
we have 
Fis iS 


Combining this result with equation (4.2), and noting that we also have H*) =o, 
we have proved the following 


Theorem 5. The affine field equations naturally admit the following system 
of equations of Maxwellian form: 


Ain =9, “Alas y=9, (4.3) 


HAGE oe Peete (4.4) 


where the Lagrangean function F is the constitutive function for this Maxwellian 
field. 

We interpret (4.2) as the constitutive equation for the exterior Maxwellian 
current S*, and we refer to equations (4.3) and (4.4) as the affine Maxwellian 
equations. 

Since A,, is a tensor, equations (4.3) insure the existence of a vector /,, such 
that 

Ia,p1 = Aap: (4.5) 


The vector 7,, so defined, will be referred to as the affine vector potential. 
Expanding equation (4.5) by means of equation (3.6), we have 


fra,p) = 22" 18, 01: 


Although the quantities L*,, are not components of a vector under general 
transformations, they do transform as components of a vector under the subgroup 
of linear transformations. Hence, under the linear group we may take L*,, as 
the affine vector potential. Under general transformations, we shall consider 
L*,, as the components of potential of the Maxwellian field tensor A,,. In 
many respects, treating L*,, as the potential for A,, is more fundamental than 
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considering a vector potential for A,,. This is due to the fact that we can find 
coordinate systems in which L*,, will vanish at a preassigned pout (7.e., a normal 
coordinate system) but not vanish for other choices of coordinate systems oe he 
freedom of choosing L‘,, to vanish at a preassigned point is equivalent to the 
freedom of choice of a reference value for the potential. 

Equation (4.4) yields, by direct computation, the topological conservation 
law for the exterior Maxwellian current 


=0. (4.6) 


If & be a scalar density, then H“* will be a contravariant tensor density, and 
the above topological conservation law becomes the tensorial conservation law 
Ste — 0; 

We have obtained results having the formal structure of MAXWELL's electro- 
dynamics by a simple separation of the effects of the tensor A,, and the tensor 
Ry» (cf. Theorem 4) without any assumptions about the form of the constitutive 
equations which describe the physical fields. This is as it should be. If we are 
to obtain representations for electromagnetic phenomena from a field theory, 
the basic field equations representative of electromagnetic phenomena (namely 
equations of Maxwellian form) should turn out to hold independently of the 
particular field which we wish to investigate; that is to say, independently of 
the particular constitutive equations used. We must, however, pick the right 
constitutive equations if theoretical predictions and experimental results are to 
coincide. Since the constitutive function for the affine Maxwellian field has 
been shown to be the Lagrangean function, we are faced with the problem which 
pervades modern physics, namely, that of finding the appropriate Lagrangean 
functions for the field equations. 


Theorem 6. If & 1s assumed a priori to be a metric space, i.e., there exists a 
tensor gg at all points of & which satisfies the postulates (2.8) through (2.10), then 
it rs impossible to obtain equations of Maxwellian form from a variational principle 
whose Lagrangean function depends only on contractions of the curvature tensor, 
the tensor gy, and its first derivatives, and the affine connection; rather, electro- 
magnetic phenomena must be introduced by the inclusion of additional field variables. 


Proof. By Theorem 3 we have that A, g=O0 since, by hypothesis, & is a 
metric space. Thus all the arguments of the Lagrangean function as hypo- 
thesized are symmetric, since C,,,;= 4A ;,—=0. There is thus no skewsymmetric 
tensor available to replace the tensor A,,, and hence additional skewsymmetric 
field variables must be introduced into the Lagrangean function in order for 
a variational principle to lead to field equations of Maxwellian form. 

The variational theory of Lanczos! might appear to contradict Theorem 6 
on first reading. This is not the case, however, if one notes that the vector Pu 
introduced by Lanczos through a “canonical transformation’? adds to the 
Lagrangean function of his theory just the skewsymmetric tensor required to 
yield equations of Maxwellian form. The result stated in Theorem 6 lies at the 


* THomas, T. Y.: op. cit., Chapter V. 


LAP AIC ZOOS Ons Electricity and General Relativity. Rev. Mod. iPhyss) 295sNome 
(July 1957). 
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heart of the problem of obtaining a field theory which naturally admits electro- 
magnetic and gravitational phenomena from a common geometric base. If we 
wish to retain the assumption of a metric, the problem of describing electromag- 
netic and gravitational phenomena from a common geometric base becomes 
formidable. The works of EINSTEIN, SCHROEDINGER and Hravat¥! presume 
an affine space having non-vanishing torsion, with its additional complexity, 
in order to represent gravitation and electromagnetism from a unified metric 
standpoint. 


5. Hamiltonian Form of the Affine Field Equations 


The study of systems of partial differential equations is significantly simplified 
in many cases if the equations are written in Hamiltonian” form. Although 
it is possible to obtain equations in Hamiltonian form by any of several methods, 
the most direct and useful for our purposes is the method of auxiliary variables 
introduced by BATEMAN and later by Lanczos}, 

The basic Lagrangean introduced in the preceding development was assumed 
to have the following functional form: 


L=L(Agg, Rag, Lay, Qo)» Ce), ¥)- (5.4) 


If we replace A,,, Rg and Q@),, by the auxiliary variables a,,, 7,3, and q.),, 
so that 
L= FL bg5,%ap, L'sy, Qe) Noy» ¥) » (5.2) 


then 4,8, %x3, L°sy, Q@), and qq, may be considered as independent provided 
the following equations of constraint are satisfied: 


Ag Hap Lop Usa, p — Cap ) 
Oe Rup <a Vx = 5 ag zie Ee) — s Spee = | Be Bee = Las al Yup» (5.4) 
0 = Gee — Qe), (5.5) 


Multiplying the equations of constraint by the Lagrangean undetermined multi- 
pliers H%*, 1% and W)® respectively, we obtain the free variational problem 


Of 2tdV —0 (5.6) 


D* 
in the variables a,,, Hee Vag fee Ly, Qe) Mee and Ww)! where 
Be (la, lng) + Oui W>_A14 (Lig ptLiep, A [Le 1%? (5.7) 


and 


KH—a ee = wba ame — De ple wis: Sle Vg) re + 


(5.8) 
ae Toe Wie’ ae L (yp, Yup Lae abs .) = 


For a brief summary see SCHROEDINGER, E.: The Final Affine Field Laws. II. 
Proc. R. I. A. 51, Sect. A. An extensive and detailed discussion can be found in 
LicHNEROWICzZ, A.: Théories Relativistes de la Gravitation et de L’Electromagnétisme. 
Paris: Masson & Cie. 1955. 

12 EpEeLeN, D. G. B.: The Invariance Group for Hamiltonian Systems of Partial 
Differential Equations I. Analysis. Arch. Rational Mech. Anal. 5 (1960). 

13 BATEMAN, H.: Phys. Rev. 38, 815—819 (1931). Lanczos, C.: op. cit. 
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Evaluating the variation of f @*dV with respect to ay, %a6» oe» respectively, 
D* . . 
and equating the results to zero in accordance with equation (5.6) yields 
ge 

He — La [*t= Loe wes — FL aye: (5.9) 
Comparing equations (5.9) with equations (3.2), (3.3) and noting that F Aap anf f 
Lop at we see that the H%* and J*°, introduced in this section as Lagran- 
Ele oii re ; ap a8 
gean undetermined multipliers, are the same as the field variables H*” and I 
introduced in the previous discussion. 


If 
det (L aug, ay) ae 0, det (L ) + 0 , det (Lae ia) yess) ai 0, 


4B Tut 


then equations (5.9) may be solved for a,,, 7g, and q{)¢ as functions of HTS, 
1*?, 1%, W)* and Qj). Using these solutions, we may eliminate a,,, 7,, and 
%o)¢ Nn equation (5.8), thus obtaining 


H= HTT LW Oey (5.10) 


Evaluating the variation of [ @*dV with respect to H%®, 1*°, L%,, W'S, 
Dt / 
Qq@), using the symmetry of L%,, I age 6L%, and the skewsymmetry of H Aes 
and equating the results to zero in accordance with equation (5.6) gives 


0= Lg. — Lap — # nb, (5.41) 
O= 3 (Lag t Lcp,0) — Lape #18, (5.12) 
0= Q@,.— # woe, (5.13) 
0=— Wel. Ho, (5.14) 
Om = 24770 26) Hs 2 0 Ore as eee ae 


Equations (5.11) through (5.15) are the Hamiltonian form of the affine field 
equations. 


Contracting equation (5.15) gives 


Satie 2H 1, at 10H". = sti 
from which we obtain 


HB” .= 7531". + 2H 74). (5.16) 
Thus, the Hamiltonian form of the exterior affine current 1s given by 
St Ol ee 2H 15). (5.17) 
Eliminating H"* , between equations (5.16) and (5.17) gives 
TM eT $0, 1 = ears oes ier noe) 


Equations (5.16) and (5.18) are equivalent to equations (5.15) and are con- 
siderably simpler to solve, since He 1K , Will be independent of H%* in general. 
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6. Conservation Laws 

The fundamental laws of physics are embodied in statements of the con- 
servation of fundamental quantities: mass, momentum, energy, charge, etc. The 
mathematical formulation of these fundamental statements involves, however, 
certain integral relations of a topological rather than a metric differential geo- 
metric nature. In the theory of general relativity and also in EINSTEIN’s unified 
field theory, conservation laws are embodied in the covariant conservation of 
a tensor 7%? of second rank, T** .=0. To obtain the fundamental topological 
conservation laws from these theories, a nontensorial quantity ¢%; must be intro- 
duced by means of the identities satisfied by the tensorial fields, which, together 
with the mixed tensor density T%, corresponding to T%®, satisfies the equation 


f (T%+t%) N,dS =o. 
This is a most important example of the topological rather than the differential 
metric structure of mathematical statements of physics. 

The present theory admits directly conservation laws of the required topo- 
logical nature. It may be proved *4 that to any system of partial differential 
equations which arise from a variational statement constructed in a Hausdorff 
space, there corresponds a collection of functions W% which satisfy the topological 
conservation laws W% ,= —0,¥. Specifically, set 


We= Lr Lane t+ Lowa 0 BL: (6.1) 

one may show that if the affine field equations are satisfied, then 
W% n= — el. (6.2) 
Using equations (3.2) and (3.3) to evaluate the terms in equation (6.1), we have 
°g = Ley, (2H™* + I™*) + QopLoe,4— Emel — OL. (6.3) 
Equations (6.2) and (6.3) constitute the basic topological conservation laws of 


this affine field theory. 


7. The Einstein Vacuum Fields as Exact Solutions 
of the Affine Field Equations in Empty Space 


We now consider the affine field equations for the case of vacuum fields. 
Definition. The space & will be said to be empty if and only if 


S*— He — wees Qe) Aus OF (7.1) 


i.e., the Maxwellian field, the exterior affine current, and the affine matter fields 
vanish at all points in &. 


For an empty &, we assume the following form for #: 
¢= constant, det (7,5) < 0, 
L(0, 74g, 0, 0, 0, 0) = 2/6 V— det (7.4) - (7.2) 
Since equation (7.2) states that # is a scalar density, the resulting field equations 


will have tensorial nature, and the colon derivative becomes the semicolon 
derivative. 


14 EpELEN, D. G. B.: op. cit., Section VII. 
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From equation (5.8) we have in this case 


aT (LH, Ug — Lh, Let tan) — 2 V— det (7,5). (7.3) 
SEE 
det (7.5) =7* (7.4) 
and 
y* 7B — Bie (7.5) 
so that rae 
LF =2C\—r*. (7.6) 
Substituting equation (7.6) into equation (5.9). gives 
jh) ae aa C7 
whence 
IP rg = A Y— 7, (7.8) 


T° 7, SAC Vr OF 
Substituting (7.8) into (7.3), we have 
H = TP (IM Lg — Lg L,4) + 2ic V—7*- (7.9) 


The only Hamiltonian affine field equations not satisfied identically under 


equation (7.1) are 
LL a pie om ea ep el el ee a 
and a o4 é€ ae 
Qf Ue OT oO Pee On (7.11) 


Contracting equation (7.11) with respect to (f, y) gives 


eee 
so that equation (7.11) is equivalent to 
eee wep (7.12) 


Differentiating equation (7.8), with respect to 1*°, we have 
2.6 = 4/0 (V—r*), 8. (7.13) 
Substituting this result into equation (7.10) and using the definition of Rag Sives 
Rag = Tap (7.14) 


in agreement with equation (5.4), as was to be expected. Substituting equation 
(7.7) into equation (7.12) gives 


(V0. (7.15) 


Equations (7.14) and (7.15) are the affine field equations for empty space. 
Expanding equation (7.15), we have 


|= Be, acai (= sy at 
or which it is sufficient to require 


V=* tuo. + top (V— 7h) ea} (7.16) 
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upon noting that Thee OF Taking determinants of VHr* 1 5=My a; we see 
that equation (7.16) implies 


A=. (7.17) 


Since 4,,—0 by hypothesis, equation (7.17) is admissible in that, by Theorem 14, 
a covariant constant scalar density field can exist in &. Thus equation (7.16) 
yields 

Teg ye Os (7.18) 


Since %,s)=0, and since det (7,,) has been assumed to be non-zero, 7, satisfies 
equations (2.8) through (2.10), and hence & reduces to a metric space under 
(7.17) and (7.18). Thus, by equation (2.15) we have 


L%3, =I, (7.19) 
since by (7.18) Q,,,=0. In this case one immediately obtains 
Ee (In = ae 


so that (7.18) implies (7.17). To this point the functions 7,, have been arbitrary 
smooth functions. Under equation (7.19), we have 


Rag = Jap (7.20) 
by equation (2.19). Substituting this result into equation (7.14), we have 
Jp =p: (7.24) 


which is a system of ten partial differential equations of the second order for 
determination of the ten potential functions 7,4. 


Theorem 7. A general class of solutions to the affine field equations im an 
empty space with Lagrangean function given by equation (7.2) 1s 


L'py= ley (7.22) 


where (1) the I%,, are Christoffel symbols based on the tensor 74g, (2) thé Yq, are 
any solution to the equations 
Seog = tap: 
Yap) = 0, (7.23) 
det (7) += 0, 


and (3) Jyg ts the third contraction of the affine curvature tensor associated with 
the connection Ig... 

By equation (2.16) 1%, is homogeneous of degree zero in the 7,, and their 
derivatives. Thus, by equations (2.5) and (2.19) J,, is homogeneous of degree 
zero in 74, and its derivatives, so that equation (7.23) is equivalent to 


[peak Diy 
B B 


é (7.24) 
Mp = Tap 
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for arbitrary 1--0 15. Equations (7.24) are just the vacuum field equations of 
EINSTEIN’s general theory of relativity with cosmological constant. Thus any 
solution to EINSTEIN’s vacuum equations with cosmological constant A gives 
an ‘7, Which when substituted into equation (7.22) constitutes an exact solution 
of the vacuum affine field equations. 

It is of interest to note that equation (7.24) exhibits the following non- 


singular solution 1%: ~~ 
—axtyAls- ava a (7-250) 


(MG Se SI = 


Thus the space resulting from the affine field equations under equation (7.1) 
is approximately Minkowskian when A is sufficiently small and x4 is restricted 
to a small enough interval. 

If we require the affine field equations to be satisfied everywhere except at 
the origin, then we obtain the solution 1” 


ty =%on = 1'e3 = — Olt = — SVG 1) aA ot 5) 
ae xv? at 42 fe Pes 


which, of course, represents empty space except in the vicinity of the origin. 
In this case, however, we obtain the three classical tests of relativity to within 
the same violations of the field equations (i.e., at the origin) as in EINSTEIN’S 
general theory of relativity. 


8. Semi-Empty Space and the Fields Resulting 
from Scalar Density Lagrangean Functions 


We now consider the affine field equations for the case of semi-empty space. 
Definition. The space & will be said to be semi-empty if and only if 


S* 20,5 hero (8.1) 


v.e., the exterior affine current and the affine matter fields vanish at all points 
of &. 

Comparing the definitions of empty space and semi-empty space, we see that 
a semi-empty space reduces to an empty space if A,,=0. 

Definition. A semi-empty space & will be said to be proper if and only if 
A,g+0 for some («, 8). We confine our attention to a proper semi-empty space. 

In the last section, we based our considerations on a particular Lagrangean 
function, partly for expediency and partly because by the very nature of empty 
space the scope of the problem was narrowed. In this section, we adopt the 
alternate course of analyzing a general class of Lagrangean functions. 


Consider the class of Lagrangean functions which are scalar densities and 
which are such that 
0), 2=0. (8.2) 


* That A must be non-zero is evident, since for 4=0 the Hamiltonian form of 


the field equations is invalid, i.e., det( RoB, Ryd) =0: 


16 McVirrir, G. C.: General Relativity and C it : 
orcieeeeee: y osmology, pp. 73. London: Chapman 


“ ToLtmMaNn, R.C.: Relativity Thermodynamics and C 1 : 
At the Clarendon Press 1934, : ST es 
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When ¥ is a member of this class we see from equations (3.2) and (3.3) that 


ae (8:3) 
and 
Hb — F sup (8.4) 


are contravariant tensor densities. The affine field equations are thus tensorial, 
and [8 = I*°,. By equations (4.2) and (8.2) we can satisfy the condition 
S* =0 only if 

I**, =0. (8.5) 
Substituting equations (8.1) and (8.5) into the affine field equations (3.7) and 


(3.8), we have 
HAS = 0, (8.6) 


IGE a0 (8.7) 


as the affine field equations in semi-empty space governed by scalar density 
Lagrangean functions. Equations (8.6) are just the second set of the Maxwell 
equations in the absence of exterior affine current. Equations (8.7) are of a 
different nature and will be seen to determine systems of potential-like functions. 


In order that & be a proper semi-empty space, we must have 
A,g=0 (8.8) 


for at least one choice of («, 8). Thus, by equation (4.5) there exists a non-zero 
affine vector potential /, == p,, such that 


fa, pe) = Aap: (8.9) 
Expanding the right side of equation (8.9) by means of equation (3.6), we have 
Tae, = 2Le(B,01- (8.10) 

Solving equation (8.10) for L’,,, we see that 
Lip =—dle+¥,p (8.11) 


is a necessary condition for & to be a proper semi-empty space. 
There are basically two ways in which equation (8.11) can be imbedded in 
a general representation for the components of affine connection: 


1) Ligy =Z'oy — x0 (fp Oy + ty 8) (8.12) 
where Z%,, is an arbitrary affinity constrained by the conditions 

Loay = Y,y3 (8.13) 
2) L%, =I%, + (fp +1, 8) + ot,8"" Sey, (8.14) 


where [;,, are Christoffel symbols based on the symmetric tensor g,g whose 
determinant does not vanish and 


Satb=—}. (8.15) 
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These alternate forms of the affinity are those which are required for represen- 
tation of the two classes of solutions to the affine field equations in semi-empty 
space. 

If 1*8 + 0, (8.6) and (8.7) are 44 equations of which only 40 are independent 
(cf. Theorem 4). Thus, if the affinity (8.12) is used, the 40 independent field 
equations together with the 4 conditions (8.13) yield 44 equations for the de- 
termination of the 44 functions Z%,, and /,. It is evident that for J%? +- 0 the use 
of (8.14) would yield an overdetermined system since there would be only 14 func- 
tions g,, and , to be determined by the 44 field equations. We have thus proved 


Theorem 8. The first class of solutions to the affine field equations in a proper 
semi-empty space with a scalar density Lagrangean function satisfying equation 
(8.2) zs given by 


1%, =L’py 10 (ts oF t h, 05) 
where the 44 functions Z%,, and },, ave any solutions to the equations 
(La, 4 ne 0 ) 
(L ryp)iy = 9, F Rap + OF 


Loup = Vy 
This class of solutions is principally of interest when we consider interacting 
electromagnetic fields and is the only possible class of solutions when the exterior 
affine current is a non-zero vector. 
The second class of solutions to equations (8.6) and (8.7) is obtained from 
the fact J*°=0 is an exact formal solution of equation (8.7). If this formal 
solution is algebraically consistent, equations (8.6) and (8.7) reduce to the 14 


equations 
A .=0, (8.16) 


Tie H mex 68 (8.17) 


Using the affinity (8.14), we obtain a deterministic system since equations (8.16) 
and (8.17) are 14 equations for the determination of the 14 unknown functions 
ae and f,. On the other hand, if (8.12) is used, we have only 14 field equations 
(8.16) and (8.17) plus the four equations (8.13) for the determination of the 
44 unknown functions 2%, and /,. The affine field equations would then only 
determine the affinity to within 30 arbitrary functions of position. It is thus 
evident that the affinity corresponding to (8.16) and (8.17) must be given by 
equation (8.14) if a deterministic system is to result. 

Since a and b in equation (8.14) are arbitrary to within equation (8.15), we 
may choose them so that equations (8.14) and (2.15) yield a unique tensor Qy4g- 
Set b=—a so that a has the value —¥# by equation (8.15). For this choice, 
equation (8.14) becomes 


L%, =I%y 3 (fp Oy 4 1, 05 1,8" Spy) - (8.18) 
Comparing equations (2.15) and (8.18), it is evident that 


Qyap aa ra bats 
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From equation (2.13) we thus obtain the following evaluation of the covariant 
derivative of g,: 


Gapsy = a ly Sap = ae Oi ep (8.19) 


and hence the field tensor R,g is expressed in terms of f,, g,g, and their deri- 
vatives by equation (2.19). The affine connection determined by equation 
(8.19), namely (8.18), is that of WrEYL’s geometry!® where g,,dx* dx" is the 
fundamental quadratic form and —+/,dx” is the fundamental linear form. 
If g,g undergoes a conformal transformation 


‘Sap = o (x) Sa p> 
then 


aay aT (Zf,0 a a) Sap — Gi a (Ino) ,) Eolas 


- Hence, if we use ’g,, in computing /%,, instead of g,, and if at the same time 
/, is transformed into 
‘Ty =f, +4 (Ino) ,, (8.20) 


we get the same affinity (8.18). This implies that the tensor g,, is fixed only 
to within an arbitrary gauge factor o(x) and that the affine vector potential /, 
is determined only to within the transformations (8.20). This arbitrariness in 
the choice of gauge is just what should be expected, since the Maxwellian field 
A,,g must be gauge invariant if it is to correspond to the field tensor of classical 
electromagnetism. We have thus proved 


Theorem 9. The second class of solutions to the affine field equations in a proper 
semi-empty space with a scalar density Lagrangean function satisfying equation 
(8.2) 1s given by 

‘py = lL, — 8 (h, 85 + hp — 8 Bey)» (8.21) 


where (1) f, ave the components of the affine vector potential, (2) 1%, are Christoffel 
symbols based on the symmetric tensor g,g, whose determinant does not vanish, 
and (3) f,, and g,pg ave any solutions to the potential equations 


H**,=0, (8.22) 
Te? = 0 (8.23) 


under which gg 1s determined to within an arbitrary gauge factor a(x) and the 
affine vector potential is determined to within the gauge transformation 


‘t, =f, +4 (Ino) ,. (8.24) 


9. Classic Vacuum Electrodynamics as Exact Solutions 
of the Affine Field Equations 


We assume & to be a proper semi-empty space and set 


La4AV Way As, he? Ww? (9.1) 


18 Wry, H.: Math. Zeitschr. 2, 384—411 (1918). SCHOUTEN: op. cit., p. 133. 
Arch. rational Mech. Anal,, Vol. 6 2 
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where h,, is a symmetric tensor of signature 2, 
b= det Vagh=9; (9.2) 
and 
WP = lea g- (9.3) 
From equations (3.2) and (3.3) we have 


He = V—h Ay, WW” (9.4) 


[Pz Oe (9.5) 


and 


Since the Lagrangean given by (9.1) is a scalar density which satisfies equation 
(8.2) and since, by equation (9.5), I*®—Q cannot yield algebraic inconsistency, 
Theorem 9 is applicable. Thus 


Ley = py — & (fy OB + te Oy — 18" Spy)» 


and the functions g,g and /, are to be determined by equations (9.4) and (9.5). 
The identical vanishing of /*? places no constraint on the choice of the functions 
Zap, and hence we may set g,g=/,, with no loss of generality. Thus, since 
the tensor h,, can be chosen at will and the tensor g,, has been shown to be the 
coefficient tensor of the fundamental quadratic form, there is no restriction on 
the choice of coordinate system for the solution of the remaining field equations 


H**,=0. (9.6) 


Equations (9.4) and (9.6) together with A,5 = ff, are just the field equations 
for the classic vacuum electromagnetic field. Moreover, the Lagrangean function 
(9.1) together with the field tensors A,, and H** are invariant under all gauge 
transformations of class C1; the Lorentz condition can be used to fix the gauge 
by solving 

vis 7, 30 (9.7) 


for (Ing) and substituting into equation (8.24); and Y hag yields the field mo- 
mentum energy tensor 

Ayn Age 2h" hePi4 — hi hr?) . 
We have thus proved 


Theorem 10. The solution manifold of the affine field equations which results 
from the Lagrangean function (9.1) in a proper semi-empty space is 1somorphic 
with the solution manifold of classic vacuum electrodynamics. 

Theorem 10 combined with the results of Section 7 shows that for one choice 
of the Lagrangean function the affine field equations describe the classic vacuum 
electromagnetic field, and for another choice they describe the classic vacuum 
gravitational field. The Lagrangean function thus plays the role of a constitutive 
function for the gravitational field as well as the electromagnetic field. This 
result is analogous to the situation in which the continuum equilibrium equations 
may describe either a solid or a fluid relative to the choice of a constitutive 
function. We thus see that the classic vacuum gravitational and electromagnetic 
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fields ave both described by the affine field equations and correspond to a simple 
separation of effects through the choice of the constitutive (Lagrangean) function 
for the particular field under investigation. 

The general problem of interacting electromagnetic and gravitational fields 
could be treated here on a formalistic basis through use of Section 8 and a study 
of various possible Lagrangean functions. An obvious first choice for # would 
be a linear combination of (7.2) and (9.1), which may be readily solved by the 
use of Theorem 9. However, this choice together with the more general forms 
of the Lagrangean function does not offer a clear cut decision as to which yields 
a more adequate picture of the true fields. It is felt that this lack of definiteness 
arises basically from the problem of a consistent and rational description of 
elementary matter. Since this problem exceeds the scope of this paper, we leave 
the examination of interacting electromagnetic and gravitational fields to a 
future communication. 
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Beitrage xu einer nichtlinearen Theorie der 5 tabilitat 
von Schichtenstrimungen langs x ylindrisch gekriimmter 
Wande gegentiber dreidimensionalen Storungen 


KLAUS KIRCHGASSNER 


Vorgelegt von H. GORTLER 


Einleitung 


Die vollstindige Erforschung der Strémungsvorgange beim laminar-turbu- 
lenten Umschlag ist ohne das Studium endlicher Stérungen nicht denkbar. 
Abgesehen von den Stérungen, deren Amplituden schon im Zeitpunkt ihrer 
Entstehung groB gegen die Grundstrémung oder vergleichbar mit ihr sind und 
daher ohne Beriicksichtigung der vollen nichtlinearen Grundgleichungen tiber- 
haupt nicht behandelt werden kénnen, gelten auch die Aussagen einer fiir kleine 
Stérungen entwickelten und bisher fast ausschlieBlich angewandten linearisierten 
Theorie im allgemeinen nur fiir kleine Zeitintervalle. Versucht man namlich, 
die Aussagen dieser Theorie auf beliebige Zeitspannen auszudehnen, so gerat 
man, wie in dieser Arbeit gezeigt werden soll, in Widerspruch zu den Ergebnissen 
einer allgemeinen Behandlung der nichtlnearen Stabilitatsphanomene. 

Erste Aussagen einer strengen nichtlinearen Stabilitatstheorie bewies J. 
SERRIN* in [/], wo in aller Strenge nachgewiesen wird, daB eine beliebige Stré- 
mung in einem beschrankten Raum gegeniiber allen denkbaren Stérungen sicher 
dann stabil ist, wenn die mit dem Maximalwert der Grundstromungsgeschwindig- 
keit und mit dem Durchmesser des raumlichen Bereichs gebildete Reynoldszahl 
kleiner als die Zahl 5,71 ist. (Hinreichendes Stabilitatskriterium.) Grund- 
str6mung und Stérungen diirfen dabei von allen drei Ortskoordinaten und der 
Zeit abhangen. Die einzigen Voraussetzungen bilden gewisse, mit den hydro- 
dynamischen Grundgleichungen stets verbundene Differenzierbarkeitseigen- 
schaften. 

Als wichtigen Anwendungsfall gab SERRIN fiir die Strémung zwischen zwei 
rotierenden koaxialen Zylindern eine im Hinblick auf die Allgemeinheit der 
zugelassenen Storungsklasse erstaunlich gute hinreichende Stabilitatsgrenze an. 
Wir kommen darauf im letzten Paragraphen eingehender zuriick. 

Eine Erweiterung des Gesichtskreises der bisherigen Stabilitatstheorie durch 
eine nichtlineare Behandlung der aufgeworfenen Probleme geben auch die in 
neuerer Zeit erschienenen theoretischen Untersuchungen von D. MEKsyn und 


* Auf die zitierte Arbeit von J. SrRRIN wurde ich dankenswerterweise von 
C, TRUESDELL hingewiesen. 
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J. T. Stuart. Sie gewinnen, nachdem sie gewisse, in ihren Konsequenzen schwer 
itbersehbare Annahmen iiber Gestalt und zeitliches Verhalten der Stérungen 
gemacht haben, durch Berechnung des nichtperiodischen Anteils der gestérten 
Grundstrémung gute Ubereinstimmung ihrer Ergebnisse mit experimentellen 
Resultaten von G. I. TaAyLor (vgl. hierzu J. T. Stuart [2]). 

Nach den Aussagen der Stabilitatstheorie kleiner Stérungen, wie sie von 
G. I. TAyLor [3], W. Torttmien [4], H. GOrTLER [5] u.a. entwickelt wurde 
(vgl. hierzu auch den zusammenfassenden Bericht von C.C. Ltn [6]), miiBten die 
St6rungen im iiberkritischen Bereich, also beim Uberschreiten eines kritischen 
Parameters (meist die Reynoldszahl) exponentiell mit der Zeit tiber alle Schranken 
wachsen. Die in den genannten Arbeiten verwendeten Stérungsansatze beschrei- 
ben jedoch das zeitliche Verhalten angefachter Stérungen nur fiir hinreichend 
kleine Zeiten richtig, wo die St6rungen noch klein gegen die Grundstrémung sind. 

Ausgehend von einer auch fiir endliche Stérungen sinnvollen Verallgemeine- 
rung der Taylor-Goértler-Wirbel werden wir, unter Beriicksichtigung der vollen 
Navier-Stokesschen Gleichungen, obere und untere Schranken fiir samtliche 
St6érungskomponenten angeben. Diese Schranken, die eine Grenze fiir das zeit- 
liche Anwachsen der St6rungen darstellen, sind endlich, falls die Grundstrémung 
selbst, die auch zustationar sein darf, fiir alle Zeiten beschrankt bleibt. 

Im letzten Paragraphen werden die hier angewandten Abschatzungs- 
methoden zur Herleitung eines hinreichenden Stabilitaétskriteriums benutzt. Das 
Ergebnis, welches allerdings nur fiir St6rungen des hier betrachteten Typs gilt, 
wird in einem Spezialfall mit einem Resultat von J. SERRIN [/] verglichen. 

Wir setzen in dieser Arbeit, neben einigen physikalisch leicht vertretbaren 
Stetigkeits- und Differenzierbarkeitseigenschaften, lediglich voraus, daB sowohl 
die Grundstr6mung wie auch die St6rungskomponenten von der Hauptstr6mungs- 
richtung unabhangig sind, wodurch insbesondere Stoérungen vom Tollmien- 
Schlichtingschen Typ ausgeschlossen werden. 

Die Untersuchungen stiitzen sich wesentlich auf den Nagumo-Westphalschen 
Satz fiir nichtlineare parabolische Differentialgleichungen, der von H. GORTLER [7] 
und K. NickeEt [8] in die Grenzschichttheorie eingefiihrt wurde. Im Anhang zu 
dieser Arbeit sind die Aussage dieses Satzes sowie einige Bezeichnungen, die sich 
bei seiner Anwendung eingebiirgert haben und hier verwendet werden, zusammen- 
gestellt. 

§1. Allgemeine Vorbemerkungen 


Wir betrachten im folgenden stets Str6mungen langs einer zylindrisch ge- 
kriimmten Flache S, die den konstanten Kriimmungsradius R besitze. Zweck- 
maBigerweise benutzen wir das in Abb. 1 dargestellte orthogonale Dreibein, 
dessen Ursprung in einen nicht genauer zu fixierenden Punkt von S gelegt wird. 
Die x-Achse zeigt dabei in Hauptstr6mungsrichtung und die y-Achse zum Kriim- 
mungsmittelpunkt hin, wahrend die z-Richtung mit einer Zylindererzeugenden 
zusammenfallt. wu, v, w bezeichnen die Geschwindigkeiten in x-, y-, z-Richtung. 

Unsere Betrachtungen beschranken sich stets auf den Bereich y2O und 
beinhalten daher zunachst nur Aussagen tiber Str6mungen langs konkav ge- 
kriimmter Wande. Die Erweiterung auf konvex gekriimmte Wande ist leicht 
durchzufiithren, interessiert jedoch in diesem Zusammenhang nicht. 


Li, Kriaus KirCHGASSNER: 


Ferner gelten unsere Uberlegungen auch im Grenzfall R=oo, also fiir die 
lings iiberstrémte ebene Platte, ohne daB wir dies jedesmal gesondert hervor- 


heben werden. ; 
Wir wollen unsere Untersuchungen auf solche Grundstromungen Uy beschran- 
ken, die Funktionen von y und ¢ allein sind, also von x und z nicht abhangen. 


(1.1) Uy = Uo (Y, t), Vp = Wy =0.- 


Fiir die Druckverteilung ergibt sich dann 


2 
(1.1) by = man! dy. 
Hierunter fallen nun insbesondere alle stationaren und instationaren Stromungen 
in der Taylorschen Versuchsanordnung, bevor die bekannte Wirbelinstabilitat 
einsetzt. Um dies einzusehen, identifi- 
ziere man einfach den auBeren Zylinder 
mit der in Abb.1 gezeichneten Flache S. 

Von der Geschwindigkeitsverteilung, 
die eine Lésung der Navier-Stokesschen 
Gleichungen ist, miissen wir fordern, 
daB sie im Innern des Strémungs- 
gebietes zweimal nach y und einmal 
nach ¢ differenzierbar ist. Die Ab- 
leitungen seien im _ abgeschlossenen 
Stromungsgebiet stetig. (Die Begriffe 
der Stetigkeit und Differenzierbarkeit 
in den stets zugrunde liegenden un- 
beschrankten Gebieten werden unten 
erklart.) 

Es erweist sich nun im folgenden als zweckmaBig, Grundgleichungen, Anfangs- 
und Randbedingungen in dimensionsloser Form zu schreiben. Daher machen 
wir von vornherein y und z mit einer fiir jedes Problem charakteristischen 
Lange h, die Geschwindigkeitskomponenten wu, v, w mit einer problemabhangigen 
Bezugsgeschwindigkeit U) und den Druck # mit @ U¢ (9 = Dichte) dimensionslos. 
Da im folgenden alle genannten Gréf8en nur in dimensionsloser Gestalt verwendet 
werden, ertibrigt sich die Einfiihrung neuer Bezeichnungen; lediglich fiir den 
dimensionslosen Kriimmungsradius setzen wir: R= R/h. 


Abb. 1. Das benutzte Koordinatensystem 


Nachdem wir angenommen haben, die Grundstrémung sei unabhaingig von 
der x- und z-Koordinate, miissen hier nur noch zwei Strémungsgebiete unter- 
schieden werden: das von zwei in z-Richtung unendlich langen, koaxialen Kreis- 
zylindern begrenzte Gebiet G,= {0< y<a; —o<z<+ oo}, (a = konst.) und das 
von der halbunendlichen ebenen Platte begrenzte Gebiet Ga {0-3 <2car 
—o<ceg<t oo}. Erwahnt sei noch, im Hinblick auf die Anwendung des Nagumo- 
Westphalschen Satzes, daB beide Gebiete sich topologisch auf das beschrankte 
Gebiet G={0<y<1; —1<z<1} abbilden lassen. 


Fiihrt man jetzt noch senkrecht zur y, z-Ebene die Zeitkoordinate ¢ ein (mit 
h|U, dimensionslos gemacht), dann bilden die Punkte (y, z,t) mit den Eigen- 
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schaften: (y, z)€G; und ¢€¢l={0<t<T} das Innere eines Zylinders Z; (i=1, 2). 
Seine Mantelflache bezeichnen wir mit S; und seine Grundflache (Schnitt mit 
der Ebene ¢=0) mit G*; S; und G¥ enthalten dabei unendlich ferne Punkte. 
Daher miiBten die Begriffe: Randwert, Stetigkeit und Differenzierbarkeit einer 
Funktion in den hier zugrunde liegenden unbeschrankten Gebieten genau definiert 
werden. Dies soll im’ Anhang dieser Arbeit geschehen, und zwar in der von 
K. NICKEL in [8] angegebenen Weise, mit der ihm die Ubertragung des Nagumo- 
Westphalschen Satzes auf nichtbeschrankte Gebiete gelang. Als wesentlich er- 
weist sich dabei nur eine Voraussetzung, da8 sich némlich das unbeschrankte 
Gebiet topologisch auf ein beschranktes Gebiet abbilden laBt. 


Nachdem dies aber fiir die hier betrachteten Gebiete gesichert ist, wie wir 
oben schon erwahnt haben, kénnen wir die oben eingefiihrten und im folgenden 
benutzten Bezeichnungen in Strenge rechtfertigen. 

Der Grundstrémung u, tiberlagern wir nun Stérungen mit den Komponenten 
Uy, Vv; und w,, die wiederum von der Hauptstrémungsrichtung unabhangig sind. 
Wir setzen also: 


u(y, 2,t) = Uy (4, t) + u(y, 2, 2) 
a4) u(y, 2, t) = V1 (¥, 2, 2) 

w(y, 2, t) = W,(¥, 2, #) 

b(y,2,b) = Po(y, t) + pr(y, 2, 4). 


Die Differentialgleichungen fiir die Storungskomponenten, die man aus (1.1) 
und den Navier-Stokesschen Gleichungen erhalt, wenn man berticksichtigt, daB 
u, und pf, diese Gleichungen befriedigen, werden erst in den folgenden Paragraphen 
angegeben, wo die Komponenten getrennt abgeschatzt werden. 

Die Anfangs- und Randbedingungen lauten, wenn u), v und w’) die Anfangs- 
verteilungen der St6rungskomponenten bezeichnen: 


w,(y, z, t) = wv (y, 2), wenn (y, 2, t) € GF, (2 ed) 
( )=0,(y, 2, t) = @,(y, 2,1) =0, wenn (VRE eS S(t 


ist. Ferner sollen die in den Differentialgleichungen auftretenden Ableitungen in 
Z, existieren und im abgeschlossenen Bereich Z,, der aus Z; durch Hinzunahme 
simtlicher Randpunkte entsteht, stetig sein. Die in §3 und §4 benutzten un- 
eigentlichen Integrale tiber G,, wo G; das durch seine Randpunkte abgeschlossene 
Gebiet G, ist, sollen gleichmaBig konvergieren, so da eine Vertauschung der 
Integrationsreihenfolge stets erlaubt ist. 


* Die folgenden Untersuchungen behalten jedoch ihre Giiltigkeit, wenn man statt 
des physikalisch sinnvollen und in (1.2) geforderten Verschwindens der Stérungen 
fiir 7-> -L oo, Periodizitat in z-Richtung verlangt, wie sie bei den linearisierten Uber- 
legungen stets auftrat. 
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§ 2. Obere und untere Schranken fur das zeitliche Anwachsen 
der ersten Stérungskomponente u, 

In diesem Paragraphen beschaftigen wir uns ausschlieBlich mit der ersten 
Stérungskomponente wu, und schatzen sie mit Hilfe des Nagumo-Westphalschen 
Satzes in einer auBerordentlich einfachen Weise nach oben und unten ab. 

Die Differentialgleichung fiir die erste Komponente 1, lautet (vgl. H. GORTLER 


[5]): 


4) Sa = 2 Llm] +( 


OUg | Uo oa | pele ) OU 
oy R-y oy k>=y 


mit dem elliptischen Operator: 


(i o? yiaih fd) 1 ee? 
~ by? Ry eOy. (Raid wale 


und mit ee (v = kinematische Zahigkeit.) Die Anfangs- und Rand- 


bedingungen wurden in den Gleichungen (1.2) angegeben. 

Wir konstruieren nun Ober- und Unterfunktionen* # bzw. u, die unabhangig 
von der Gestalt der Stérungskomponenten v, und w, sind. Zu diesem Zweck 
miissen wir von # und # verlangen, daB sie Lésungen der Differentialgleichungen 
Ou uU OUp Uo ou 

= —— — ° —- 0) 
Oy R—-y oy a Ry Oz 


(2.2) 


sind, wodurch das Verschwinden der Koeffizienten von v, und w, erreicht wird. 


Die allgemeine Losung der Gleichungen (2.2) lautet: 


U(y,t) = — u(y, t) + —— 
(2.3) 
%(y,t) = — u(y, t) + 


e 
he 


Wir kénnen nun iiber « und « so verfiigen, daB #@ und % die folgenden, eine Ober- 
und Unterfunktion kennzeichnenden Ungleichungen befriedigen: 


Ou 1 


vesnalcs 
(2.44) ee ans 
U 1 ~ 
Gra RAE 
und 
ti(y,t) > w(y,z), falls (y, 2, t)€G#; 
(2.4b) ae t) > 0, falls (y,z,t)€S,; 
a(y,t)<u(y,z), falls (y,z, t) CGF; 
u(y, t) <0, falls (y,z,t)€S,, (= 1,2). 


Falls die betrachtete Stérung periodisch in 2z-Richtung ist, braucht man die 
Ungleichungen nur auf G* und fiir ein Periodenintervall zu fordern. 


ro dha Zusammenhang mit der Entwicklung des Nagumo-Westphalschen Satzes 
und verwandter Satze haben sich die Bezeichnungen ,,Ober-“ und ,,Unterfunktion“ 
eingebiirgert. Ihre genaue Definition wird im Anhang wiedergegeben. 
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Fiir « und & setzen wir an: 


(2.5) 


~ 


Fiihrt man diese Beziehungen in (2.3) und # bzw. @ in (2.4a) und (2.4b) ein, 
so zeigt sich, daB die genannten Ungleichungen erfiillt sind, wenn man 


Reese he ieltedllne 
UR — 9) Moly ones tei € 
a Wh aoe «| 


(y, 2, ¢) 


a* — Max | 
(2.6) 


(1%) [[(R — y) uo(y, Dlg Wess ao eee 
é>0O, sonst beliebig und B>0: p> 0, sonst beliebig wahlt. 
Aus dem Nagumo-Westphalschen Satz (vgl. Anhang) folgt dann: 
ar ob! a ob! 
— U(¥,t) + Rady, < Uy (Y,2,t) << — Uq(y, t) + fey . 


LaBt man nun in diesen Ungleichungen B und B und in (2.6) e gegen Null 
gehen, so folgt: 


(2.7) ~ (9) +a 


(od 


Rk 


S u(y, 2,4) S — u(y, 8) + ieee 


Y) 
mit 


[(R — ¥) u(y, lo,nes 0 
{i a AhleN Cr agree ool) i 


ee eat "| 


[(R — y) u(y, t) ly,€ Se 0) 

Die Beziehung (2.7) stellt nun die angekiindigte Abschatzung der ersten 
Stérungskomponente wu, dar. Sie gilt, weil wir tiber 7 keine einschrankenden 
Annahmen gemacht haben, fiir jedes zeitliche Intervall OS¢<7, und damit 
fiir alle Zeiten ¢: OS¢t<oco. Dabei forderten wir von u, lediglich gewisse Stetig- 
keits- und Differenzierbarkeitseigenschaften und ferner, daB es im Unendlichen 
verschwindet oder in z-Richtung periodisch ist. 

Ist die St6rungskomponente u, zum Zeitpunkt ¢=0 ihrer Erzeugung klein 
gegeniiber der Grundstrémung u(y, [=0), so hegen die hier abgeleiteten oberen 
und unteren Schranken von uw, in der GroBenordnung von u(y, ¢). Ist dagegen 
u(y, 2) vergleichbar mit u(y, f=0), so setzen sich die obere und die untere 
Schranke aus dem Maximal- bzw. Minimalwert der Grundstr6mung und den 
entsprechenden Werten der Anfangsverteilung u(y, z) zusammen. Wenn u(y, @) 
selbst fiir alle Zeiten beschrankt ist, was insbesondere fiir alle stationaren Grund- 
stromungen zutrifft, gilt dies nach (2.7) auch fiir die St6rungskomponente %,. 

Zum Schlusse soll die Scharfe der Abschatzung (2.7) an Hand eines wichtigen 
Beispiels gepriift werden. 


,t) 


* Um dieselbe Abschatzung fiir den Fall R=oo durchfiihren zu kénnen, hat man 
nur @/(R—y) bzw. %/(R—y) durch zwei Konstanten ¢ bzw. ¢ zu ersetzen. 
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Up sei die stationare ebene Strémung in der Taylorschen Versuchsanordnung ; 
@, und @, seien die Winkelgeschwindigkeiten des inneren und des duBeren 
Zylinders, R, und R, ihre Kriimmungsradien. Als charakteristische Lange / 
wahlen wir die Spaltbreite h=R,—k,, als Bezugsgeschwindigkeit die Umdre- 
hungsgeschwindigkeit des inneren Zylinders Uy= Ry @. 


Setzt man nun RoR Re R—1=R,h=R,, 


so hat uw, die folgende Gestalt: 


B 
mit i i? — (©, /a,) Fes = Fa R, R3(1 — w/o) ) 
Rite Rg 2R,—1 


Um die Konstanten & und & in den Gleichungen (2.7) zu berechnen, haben 
wir den in den geschweiften Klammern stehenden Ausdruck nach oben und 
unten abzuschatzen. Zu diesem Zweck betrachten wir die Funktion 


Ply) = (Ry — ¥) Uo (¥) - 
Sie wachst fiir 4 <0 monoton im Intervall 0<y<1, fallt monoton fir 4>0 
in diesem Intervall und ist fiir d=0 konstant. Daher gilt: 
_ [Rup +R. (1), tire. A <0; 
R, ue + Rau, (0), tae VAs 0: 
a—R,u®; A beliebig, 


mit 
uy = Max u(y, z) 
Gs 


us) = Min u(y, z). 
Gi 


Wir schatzen nun noch die rechten Seiten der ersten Beziehung (2.7) nach 
oben ab. Fir A <0 folgt: B>R,. Beachtet man noch, daB m#(41)=4 und (0) = 


Ry W,/R, w, ist, so zeigt man leicht, daB die Funktion 
= ee i 
P1(Y) u(y) 4 Fils 
im Intervall OX y<1 monoton fallt, so daB sich hieraus, zusammen mit (2.7) 
die Abschatzung ergibt: 


(2.8) — to(y) +t ul) < aa (y, 2,1) S— 


) 


(0) 1 R2 R2 
uU ——— Rio ROSS 
pease Bs j eae ere: ai 202] 


fiir A<0, d.h. w,/o,< R?/R3. 


Entsprechend kann fiir A>0 verfahren werden. In diesem Fall wachst 
namlich die Funktion 


P2(¥) = — uo(y) +g? — 149 (0) 
R.—-y 
monoton im Intervall O< y<1. Daher gilt: 


R = 
(2a = tee) pt yl) aay (yy gp) ta 0) ee a 
0 R,—y = L(V ) ) = R ay Uy aP Fes Or [R53 @» ky 4] 


fiir A>0, dh. w/a, > R2/R2. 
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Aus den oberen Schranken fiir u, in den Beziehungen (2.8) und (2.8a) erkennt 
man, daB diese Schranken immer kleiner werden, wenn A gegen Null geht. 
Bei gegebener Rotationsgeschwindigkeit des inneren Zylinders (gegebene Reynolds- 
zahl in (2.1)) und geeigneter Wahl der Rotationsgeschwindigkeit des auBeren 
Zylinders (A hinreichend klein), kann man dann erreichen, daB sich die obere 


Schranke fiir u, um beliebig wenig von dem Term R,/(R,— y) uv) unterscheidet. 


Endliche Stérungskomponenten sind damit fiir geniigend kleine A-Werte im 
wesentlichen durch den Maximalwert ihrer Anfangsverteilung nach oben be- 
schrankt. 

Fiir kleine St6rungskomponenten lassen sich aber aus dieser Beziehung noch 
weitere Schliisse ziehen. Da namlich im Rahmen der linearisierten Stabilitats- 
theorie, die ja das Verhalten hinreichend kleiner Stérungen angendhert richtig 
beschreibt, das Vorzeichen von u, nicht festgelegt ist (uw, ist Lésung eines homo- 
genen Randwertproblems), stellt die rechte Seite von (2.8) und (2.8a) sogar eine 
obere Schranke fiir |u,| dar. Liegt «°}, nun in dem durch die lineare Theorie 
erfaBbaren Bereich, so lehrt die Beziehung (2.8a), daB |u,| fiir geniigend kleine 
A-Werte niemals aus diesem Bereich heraustreten kann. u, ist demnach sicher 
nicht angefacht. Die Grundstrémung wird also gegen kleine St6rungen immer 
stabiler, die kritische Reynoldszahl immer gréBer, wenn man sich dem Wert 
A=0 nahert. Diese Aussage, die schon im Rayleighschen Kriterium enthalten 
ist, ergibt sich hier direkt aus den Abschatzungen. 

Die vollstandige Aussage dieses Kriteriums, da8 namlich die Grundstr6mung u, 
_ sogar fiir AO gegen kleine Stdrungen stabil ist, laBt sich jedoch aus den 
obigen Abschatzungen nicht beweisen. 


§ 3. Schranken fiir das zeitliche Anwachsen der St6rungskomponentenr, undw, 

Da die Abschatzung der Stérungskomponenten v, und w, mit Hilfe des im 
vorigen Paragraphen verwendeten Nagumo-Westphalschen Satzes nur grobe 
Schranken liefert — man kann lediglich zeigen, daB v, und w, héchstens linear 
— mit der Zeit anwachsen —, werden wir in diesem Paragraphen zu einer anderen 
Methode greifen. Wir gehen von den Grundgleichungen iiber zu den Gleichungen 
fiir die gesamte, im durchstrémten Raum vorhandene kinetische Energie der 
Komponenten 7, und w,. Auf dem Ergebnis des vorigen Paragraphen aufbauend, 
kann man dann mit einer schwachen und physikalisch belanglosen Einschrankung 
zeigen, daB dieser Energiebetrag fiir alle Zeiten unterhalb einer festen, von der 
Zeit unabhangigen Schranke bleibt, wenn nur die Grundstrémung selbst be- 
schrankt ist. Dies bedeutet aber die Beschranktheit fiir die Komponenten selbst. 

Die Differentialgleichungen fiir die betrachteten St6rungskomponenten erhalt 
man wiederum durch Einsetzen von (1.1) in die Navier-Stokesschen Gleichungen. 
Wenn man beachtet, daB u, und ~, diese Gleichungen schon befriedigen, folgt 
dann: 


Oy Milos g ae Op, Uy (2Uy +) OY a Ov, v1, 
t Re = ey R—y ey Oz 
OW, 1 Coe Op, Ow, . OW, 
ff v,| => — — — 4, — —— 
(3-1) ot Re L,[“1) OZ oy : 1 ENG 
OU; 0, OW, 
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Multipliziert man nun die erste Gleichung in (3.1) mit (R—y)v, und die zweite 
Gleichung mit (R —y)w,, integriert man dann tiber den durchstromten Raum 
und addiert man schlieBlich die erhaltenen Ausdriicke, so ergibt sich, nach 
Durchfiihrung mehrerer partieller Integrationen und unter Beriicksichtigung der 
Kontinuitatsgleichung : 


5 yf R—¥) ott wh dy det 


2 bt 
G; Lat, ; 
62) 46. fR—) [(h + (Bl + ee + (GS) + er 
g, 
ee [0 14 (2t4g + m4) dy dz, 
G, 


dabei bedeutet G, das durch seine Randpunkte abgeschlossene Gebiet G,. Setzt 
man voriibergehend: v,= 7, sin z und w,=@, sin z, so folgt mit Hilfe dieser Sub- 
stitution: 
(3.3) {, (R — y) (Sy) ay dz =[(R — y)vidydz + f(y sin? zdydz 

G, G, G, 
und eine entsprechende Beziehung fiir w, und @,. 

Wenn man nun (3.3) in (3.2) einsetzt, die rechte Seite von (3.3) mit der 
Cauchy-Schwarzschen Ungleichung nach oben und die linke Seite durch Ver- 
nachlassigung der Quadrate der Ableitungen sowie des Terms v/(R — yy)? nach 
unten abschatzt, so folgt: 


5 | RN Wit wh dydz+ sf (R—y) (i+ wh) dy dz 
(3.4) iz 


Gy 
<| [eidyde. [iam +1m)%dy dal, 


G G 
Fiihrt man noch die Hilfsgr6Ben 


[pQP=J(R—y) i+ widydz; [pO P=fuk(2uytm)*dydz 
G, G, 
ein, so erhalt man aus (3.4) die einfache Beziehung: 


dit 1 
(3.5) EE a zai 
Aus ihr laBt sich nun eine obere Schranke fiir y ableiten, wenn nur y fir alle 
Zeiten beschrankt ist. Hierfiir miissen wir aber zunichst, wie im vorigen Para- 
graphen, voraussetzen, daB a» fiir {—>co endlich bleibt. Ferner muB gefordert 
werden, daB das Integral: 
(3.6) [aidy dz, 


Gy 


vs Falls die Storungskomponenten periodisch in z-Richtung sind, umfa8t der 
Integrationsbereich nur eine Periode. 
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mi - 
: U,(¥, 2) = Max ilma(y, z,)]} 
existiert. Nachdem wir aber im vorigen Paragraphen gezeigt haben, daB |u| 
in jedem Punkt (y, 2, 4) €Z, beschrankt ist, erscheint diese Annahme vom physi- 
kalischen Standpunkt aus nicht einschneidend. Betrachtet man etwa den Taylor- 
schen Fall der koaxialen Zylinder, so ist die obige Forderung schon dann erfiillt, 
wenn man von der idealisierten Vorstellung unendlich langer Zylinder abgeht 
und endliche Zylinderhéhen annimmt. 
Nach den Ungleichungen (2.7) 14Bt sich | w)-+-u,| durch Max {4/(R — y) (&, |&|)} 
nach oben abschatzen, woraus sich, zusammen mit der Annahme (3.6), die 
folgende Ungleichung ergibt: 


(3:7) |y@)| Ss | fatay az) Max {|u| + |u% +u%,|}< M; (M = konst.). 
es 1 


Nun folgt aber aus (3.5) und (3.7): 
(3.8) p(t) < Nel-VRet + Re M, 


wobei 
— 5 


NZ{ f(R—y) [0 + (w)] dy az} — Rem 


Gy 


ist. Die Beziehung (3.8) ist aber unsere Behauptung, daB namlich die Komponenten 
aller betrachteten Stérungen fiir alle Zeiten unterhalb einer festen Schranke 
bleiben *. : 

Um analoge Abschatzungen fiir Roo zu erhalten, hat man die Gleichung 
(3.4). fir R=co mit v, und w, zu multiplizieren, iiber G, zu integrieren und die 
so erhaltenen Beziehungen zu addieren. Die Beschranktheit der Stérungs- 
komponenten v, und w, ergibt sich dann, wie dies oben im Falle der zylindrisch 
gekriimmten Wand dargelegt wurde. 


§ 4. Ein allgemeines Stabilitatskriterium 

Die im vorigen Paragraphen benutzte Methode soll zum SchluB zur Her- 
leitung des in der Einleitung angekiindigten Stabilitatskriteriums verwendet 
werden. Wir schranken dabei die Allgemeinheit der Aussage wiederum nur durch 
die Forderung ein, daB die Grundstrémung und samtliche Stérungskomponenten 
von der Hauptstromungsrichtung unabhangig sind. 

Die Beziehung (2.1) wird jetzt noch mit (R—y) u, multipliziert, sodann 
iiber den Bereich G, integriert, und endlich wird die so erhaltene Gleichung zur 
Gleichung (3.2) addiert. Benutzt man wiederum (3.3) und vernachlassigt danach 
alle Quadrate der Ableitungen, so ergibt sich mit Hilfe der Cauchy-Schwarzschen 
Ungleichung: 


(4.1) 1 IPO + 2 POPS CIPO? 


* Yom rein mathematischen Standpunkt aus mtiBte diese Aussage etwas modi- 
fiziert werden. Aus der Beschranktheit der Gesamtenergie y folgt namlich lediglich, 
daB die Funktionen v, und w, ,,fast iiberall‘‘ beschrankt sind. 


30 Kraus KrrRCHGASSNER: 


it (pOP=S (Ry) e+ 2 + w}) dy de 
G 
und 
ah ried eae 
Os a Max { oy i R—y i 


Aus (4.4) folgt nun aber sofort, daB fiir Re<1/C 

lim |G (¢)|? =0 

t— oo 
gilt. Mithin sind alle infinitesimalen und alle endlichen Stérungen im Bereich 
Re<1/C gedampft, wenn nur die Grundstroémung und die Storungen von der 
Hauptstromungsrichtung unabhangig sind. 

Eine einfache Rechnung ergibt im Taylor-Fall: 
Se 
Ry 

Alle stationaren Stromungen in der Taylorschen Versuchsanordnung sind daher 
im Bereich 


R, (R3—R3) 1 { 3 1 ! 
5 R IO a ea NS —— ph ee (a Aeeels 
‘ a) ee R3|1—a@,/@,| |1—a,/@, | R,/h . (F2/h)? 


gegentiber von x unabhangigen Stérungen beliebiger GréBe stabil. (Die Reynolds- 
zahl Re wird hierbei wie in §2 mit der Rotationsgeschwindigkeit des inneren 
Zylinders R,@, und mit der Spaltbreite h=R,—R, gebildet; Re=R, a, hy.) 
Ein Vergleich der in (4.2) abgeleiteten Stabilitatsgrenze Re, mit der von 
J. SERRIN in [7] angegebenen (wir nennen sie Fe,) fallt zugunsten der letzteren 
aus. SERRINS Ergebnis lautet namlich, in den hier verwendeten Symbolen ge- 
schrieben: 
ae a (R3— Ri) 
R, R3 (log (R,| R,)]? |1—@2/a, | 


Re, 


Es gilt stets 2* Rey<Re,, wobei das Gleichheitszeichen im Grenzfall R,/h—>oo 
angenommen wird. Da SERRINs Stabilitatsgrenze dariiber hinaus auch fiir x- 
abhangige Stérungen gilt, ist die Aussage der Beziehung (4.2) schon vdllig im 
Serrinschen Ergebnis enthalten. 

Als zweites Beispiel werde die durch ein Druckgefalle erzeugte stationdre 
Schichtenstrémung zwischen zwei koaxialen Zylinderwianden untersucht. Die 
Spaltbreite h=R,—R, sei wiederum als charakteristische Lange, die Maximal- 
geschwindigkeit der Grundstrémung uy als Bezugsgeschwindigkeit gewahlt. u, 
hat dann die Gestalt: 


a b end = 
Uo(¥) =y |a(Ry— y) 4 pee (R, y) log(R,— y)|, 
mit es oy als od 
qa — Plog R,—Rj log ky 
RR 
a ee, Rs 
p—— FAM og Be 
Ria oe 
Nas Ry 
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dabei ist Ry=R.— Yo, und yp ist durch die Gleichung up (V9) =, bestimmt. 
Der oben abgeleitete Stabilitatsbereich Re<1/C lautet in diesem Fall: 


2RY|2b+ RG] 


AR: R Re S== = 
>) ea iad AD aT 


Da sich 2b/R3 fiir groBes R, wie 2b/R? verhalt, ergibt sich an Stelle des Kriteriums 
(4.3) die einfachere Beziehung: 


(4.3.a) Re <2R,. 


Der Stabilitatsbereich der betrachteten Kanalstrémung wachst also fiir hin- 
reichend groBes R, monoton mit wachsendem Ra Die ebene Kanalstromung 
selbst ist véllig stabil gegeniiber Storungen des hier betrachteten Typs. 

Ein Vergleich mit der von J. SERRIN fiir beliebige Grundstrémungen ange- 
gebenen Stabilitatsgrenze Re=5,71 — sie wird in der zitierten Arbeit allerdings 
nur fiir beschrankte Bereiche bewiesen — zeigt, da die Beziehungen (4.3) und 
(4.3a) fiir geniigend groBes R, gréBere Stabilitatsbereiche liefern. Dies war zu 
erwarten; denn die ebene Kanalstrémung wird bekanntlich schon gegeniiber einer 
speziellen St6rungsart der bei SERRIN zugelassenen allgemeinen St6rungen, nim- 
lich der Tollmien-Schlichting-Wellen, bei endlichen Reynoldszahlen instabil. 

Im Falle der ebenen Platte (Roo) lautet das Stabilitatskriterium entspre- 


chend: Fiir Re<1/C, mit C=Max {4 Ee 
Ze BN ODM 


endlichen, von x unabhangigen Stdrungen gedampft. 


} sind alle infinitesimalen und alle 


§5. Anhang 


Der Nagumo-Westphalsche Satz. Wir haben in der vorliegenden Arbeit den 
Nagumo-Westphalschen Satz an einer wesentlichen Stelle, namlich bei der Ab- 
schatzung der St6rungskomponente wu, benutzt. Desgleichen verwendeten wir in 
diesem Zusammenhang die Bezeichnungen ,,Ober“‘- und ,,Unterfunktion“‘. Diese 
Bezeichnungen sowie die Aussage des Satzes sollen hier im Anhang kurz wieder- 
gegeben werden. Bei der Formulierung der Voraussetzungen beschranken wir 
uns auf die in dieser Arbeit behandelten speziellen Str6mungsgebiete und Rand- 
bedingungen und verweisen im itibrigen, was den Beweis und die Erweiterung 
auf allgemeinere Bereiche anbetrifft, auf die Arbeit von K. NICKEL [8]. 

Im zweidimensionalen euklidischen Raum, in dem wir die kartesischen 
Koordinaten y, z eingefiihrt denken, sei ein Gebiet G={0< y<1; —1<z<+1} 
gegeben; ferner ein Intervall J={0<t<T}. Die Gesamtheit der Punkte (y, z, ?) 
mit (y, z)€G und ¢€J bilden das Innere eines Zylinders Z, dessen Mantelflache 
wir mit S und dessen Grundflaéche (Schnitt mit der Ebene ¢=0) wir mit G* be- 
zeichnen. Z sei derjenige abgeschlossene Bereich, der aus Z durch Hinzunahme 
samtlicher Randpunkte hervorgeht. Ferner sei die Funktion u(y, z,¢) Lésung 
des folgenden Randwertproblems: 


Ou du ou, Ou = 
“Ot OV) MC Lian CN ee 
(5.4) u(y, Z,¢=0) = u(y, z); 
u(y, Z,t) =0, falls (¥,2,2)ES 


=F (i, Vi 2 W} 
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uw existiere und sei stetig in Z. Die ersten und zweiten Ableitungen von uw nach 
y und z und die erste Ableitung von u nach ¢ sollen in Z existieren und inZUG*US 
stetig sein. oe 

Die Funktion F sei eindeutig definiert und beschrankt fiir (y, z, 4)€Z und 
fiir beliebige reelle Werte der restlichen Variablen. Dariiber hinaus soll F in den 
zwei letzten Variablen schwach monoton anwachsen. Es sei also 


F(t, 2,03 > Qa; Qt %01) Sb Ys 23 Gi, G23 Na» Ia2) * 


falls %1S 42 und Jo1S Joe gilt. 

ai bzw. wv heiBen Ober- bzw. Unterfunktion von uw, wenn sie allen an u gestellten 
Stetigkeits- und Differenzierbarkeitsforderungen gentigen und die folgenden Un- 
gleichungen erfiillen: 


ge en 
(5.1) i (y, Z,¢ = 0) > u (y, 2); 

“(y,2,t)>0, falls (¥,2;2);E 5. 
bzw. 

ee 
(5.1 b) “(Y, 2, t= 0) < wO(yr2)5 


u(y,z,t)<0, falls (y,2,0€S. 
Es gilt dann der 
Satz von Nagumo-Westphal. Ist die Funktion w eine Lésung des Randwert- 
problems (5.1), die sdmtliche oben angefiihrten Stetigkeits- und Differenzier- 
barkeitseigenschaften besitzt, und ist # eine Oberfunktion und # eine Unter- 
funktion von w, dann gilt im ganzen Béreich Z: 


ty, 2,t) < u(y, 2,t) <<u(y,z,t). 


Die Erweiterung dieses Satzes auf nichtbeschrankte Gebiete kann nach 
K. NrcKE1 [8] in einfacher Weise geschehen. Wir fithren den Gedankengang an 
Hand der in §1 definierten unbeschrankten Gebiete G; durch. Zundchst werden 
diese Gebiete topologisch auf G abgebildet. Dies leisten die Transformationen: 


* 
= ayes a (=A eet); fire 
Ss ae 
sdemar pear Miia caaer iy Wn eS UE G7 


Die Randpunkte von G; gehen dabei in Randpunkte von G iiber; daneben enthilt 
aber der Rand von G noch weitere Punkte (etwa z*=-+-1), die als Bildpunkte 
der unendlich fernen Punkte von G; aufgefaBt werden. 

Durch die Transformationen sind nun auch fiir das Gebiet G; und damit 
auch fiir Z; alle Randpunkte durch ihre topologischen Bilder eindeutig erklart. 
Definiert man nun noch ,,Randwerte‘‘ einer Funktion in den unendlich fernen 
Punkten durch die Werte der Funktion in den zugeordneten Punkten, ebenso die 


* Man erkennt leicht, daB die rechte Seite der Gleichung (2.1) die hier an F 
gestellten Forderungen erfiillt. 
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Begriffe: Konvergenz einer Punktfolge, Stetigkeit und Differenzierbarkeit einer 
Funktion im Bereich G; durch die entsprechenden Eigenschaften in G, dann 
laBt sich die Aussage des Nagumo-Westphalschen Satzes auch fiir die hier ver- 
wendeten nichtbeschrankten Grundgebiete beweisen. 


Zusammenfassung 


In dieser Arbeit wird das Verhalten dreidimensionaler Stérungen einer statio- 
naren oder instationaren Schichtenstrémung langs einer zylindrisch gekriimmten 
Wand untersucht. Grundstro6mung und Stérungen sollen von der Hauptstr6- 
mungsrichtung unabhangig sein, unterliegen aber sonst keinerlei physikalisch 
einschrankenden Voraussetzungen; insbesondere werden Stérungen_beliebiger 
Gr6éBe zugelassen, wobei die bekannten Taylor-Gortler-Wirbel als Spezialfall in 
diesen Untersuchungen mit enthalten sind. Es werden Schranken angegeben, 
liber die hinaus die Stérungskomponenten niemals anwachsen kénnen._ Diese 
Schranken bleiben fiir alle Zeiten endlich, wenn die Grundstrémungsgeschwindig- 
keit selbst endlich bleibt. 

Ferner wird gezeigt, daB alle betrachteten infinitesimaien oder endlichen 
Stérungen gedampft sind, falls die Reynoldszahl einen gewissen, wiederum von 
der Grundstr6mung abhangigen und oben angegebenen Wert nicht iiberschreitet. 

Diese Abhandlung aus dem Institut fiir Angewandte Mathematik der Universitat 
Freiburg 1. Br. und dem Institut fiir Angewandte Mathematik und Mechanik der 
Deutschen Versuchsanstalt fiir Luftfahrt, Freiburg 1. Br., ist der 2. Teil meiner 
Dissertation: Beitrage zur hydrodynamischen Stabilitatstheorie. Freiburg i. Br., Mai 


1959. Die Arbeit wurde von dem Wirtschaftsministerium des Landes Baden-Wiirttem- 
berg gefordert. 
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1. Introduction 


The present paper is concerned with the theory of integration of the displace- 
ment-equations of motion in the linear theory of homogeneous and isotropic 
elastic solids. These equations, which in their general form originated with 
Caucuy [1]! (1828), in vector notation, and for the case of vanishing body forces, 
appear as 


1 _@ Ou 
Vat ie Soe ou ay Wane (4) 
Here u(x, ¢) is the displacement vector? — x standing for the triplet of rectangular 
cartesian coordinates (1, %2, X¥3) — ¢ denotes the time, the constants @, mw, and » 
designate the mass density, the shear modulus, and Polsson’s ratio, respectively, 
while V is the usual del-operator. 

In his mémoire on the equilibrium and motion of elastic bodies, PoIsson [2] 
(1829) introduced a class of particular solutions to (1) that is obtained by assuming 
the displacement vector to be the gradient of a scalar potential. In a subsequent 
addition to the foregoing paper, Portsson [3] (1829) proved? what appears to 
be the earliest theorem concerning the general integral of the equations of motion. 
He showed that every (sufficiently regular) solution of (1) admits the represen- 
tation 


u(x,t) =u'(x,t)+u"(x,t), w=Vo, V-u’=0, (2) 
where y(x,?) and w’’(x, t) satisfy the wave equations 


ip=0, ([u"=0, (3) 


* Numbers in brackets refer to the list of references at the end of the paper. 

* Throughout this paper letters in boldface represent vectors. The symbols ‘:”’ 
and ‘‘ x * are used to indicate scalar and vector multiplication, respectively. 

* This proof is summarized by TopHUNTER & PEARSON [4], Vol. 3p. 273 et seg: 
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in whic A=? Be (n = 1, 2), 
FA 4 
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Potsson thus established the fundamental fact that the complete solution of (1) 
may be expressed as a superposition of an irrotational and an equivoluminal 
motion, associated with the respective wave velocities c, and cy. The results 
contained in [2], [3] rest on the molecular hypothesis of elastic action which, 
in the case of the isotropic medium, implies the existence of but a single inde- 
pendent elastic constant and is equivalent to choosing y= +4 in (1). Consequently, 
the wave velocities c,, cz defined in (4) were found by Poisson to obey the relation 
cj=3c3. His proof [3] of the preceding theorem, however, in no way depends 
on this particular choice of the parameter ». 

STOKES [5] (1851), in dealing with the general isotropic medium of unlimited 
extent, approached the integration of the equations of motion by deducing from 
(1) that the dilatation # and the rotation w, which are given by 


Ot ie Ve yx, 2) == VX Ue (5) 
must satisfy the wave equations 
ie 0 le 0 (6) 


He then made use of an integral representation for the displacement u(x, ¢) in 
terms of # and w. 

The general solution (2), (3), due to Poisson, does not involve the vector 
potential appropriate to the solenoidal displacement component wu”. A solution 
of (1) which employs both a scalar and a vector potential apparently made its 
first appearance‘ in LAME’s [6] (1852) treatise. LAME noted that every vector 
field of the form 


u(x,t) =Vpep+tVxp (7) 
meets (1) provided (x, ¢) and (x, ¢) are particular solutions of 
Hig=0, Lib, (8) 


as is readily confirmed by substitution. 


The question as to the completeness of the solution (7), (8) was raised by 
CLEBSCH [7] (1863), who asserted that every solution of the equations of motion 
(1) in turn admits the representation (7), (8). With a view toward proving this 
theorem, CLEBSCH independently demonstrated the completeness of the solution 
(2), (3) by an argument which is somewhat less direct than that originally pro- 
duced by Potsson [3], but without specializing the value of y; he subsequently 
assumed without proof that every w’’(x, 7?) satisfying the last of (2) and the 
second of (3) may be represented by 


w(x) =Txp, O3p=o. (0) 
Another inconclusive proof of CLEBSCH’s theorem was given later by KELvIn [8] 
(1884). 


4 See TODHUNTER & PEARSON [4], Vol. 2, Pt. 2, p. 99. 
5 See Lecture IV, p. 41 et seg. This lecture was delivered in 1884 although [8] 
was published in 1904. 
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Dune [9] (1898), after expressing reservations about CLEBSCH'S [7] com- 
pleteness proof®, supplied a rigorous proof of the theorem previously announced 
in [7]. In a second paper, DuHEM [10] (1900) extended the result established 
in [9] to more general systems of linear partial differential equations. . 

Unfortunately, DunEm’s [9] original paper, which is not readily accessible, 
appears to be little known, and treatises on the theory of elasticity usually 
refrain from discussing the generality of Lams&’s solution (7), (8), or take its 
completeness for granted. In contrast, SNEDDON & BERRY [127]? (1958) state 
that Lamé’s solution is incomplete if the region occupied by the medium has a 
boundary. A proof of CLEBSCH’s theorem appearing in PEARson’s [12] (1959) 
book is open to objections, as is a brief argument advanced in this connection 
by BisHop [13] (1953). 

In view of the preceding remarks, and because of the far-reaching importance 
of Lam&’s solution in elastokinetics, it seems appropriate to include DUHEM’s 
proof of CrEBscH’s theorem in this paper (Section 2). At the same time, the 
present vectorial version of the proof is apt to bring out more clearly the essential 
simplicity of DUHEM’s reasoning. 

In Section 3 we explore the significance of LAME’s solution within the equi- 
librium theory. While there is nothing in DUHEM’s completeness proof that 
impairs its validity when u happens to be independent of the time, (7) and (8), 
with y and & functions of position alone, clearly cannot be the general solution 
of (1) in the equilibrium case: for in this event (8) would degenerate into Laplace 
equations and hence, by (7), every equilibrium solution u(x) of (1) would have 
to be harmonic (rather than merely biharmonic). This apparent contradiction 
is resolved by observing that even in the equilibrium case the potentials m and, 
in general, remain time-dependent. Thus every equilibrium field u(x) conforming 
to (4) may be represented as the sum of an irrotational and an equivoluminal 
motion, the time-dependence of which cancels upon addition. In the section 
here referred to we determine the explicit manner in which (x, t) and (x, 2) 
involve ¢ when u= u(x) by establishing a connection between the LAME potentials 
y, ~ and the PAPKOVICH-NEUBER [/4], [15] potentials of the equilibrium theory. 
The resolution of an equilibrium solution into irrotational and equivoluminal 
dynamic components is illustrated with a specific example. This type of resolution 
has an analogue in the theory of the repeated wave equation, as is made clear 
at the end of Section 3. 

Section 4 is concerned primarily with the relation between Lamé&’s solution 
(7), (8) and a solution of (1) discovered by IACovACHE [16] on the basis of a formal 
operational scheme due to G. C. Motstt [17]. We show that IACOVACHE’S solution, 
which is the dynamic counterpart of GALERKIN’s [78] general equilibrium solution, 
may be reduced to Lamg&’s solution with the aid of a simple transformation. 
Since the completeness of Iacovacne’s result was established in a previous 
paper [19], this transformation supplies an alternative proof of CLEBSCH’s com- 
pleteness theorem for Lamx&’s solution. At the end of Section 4, a recent general 
solution due to NoLt [20] of the equations of motion for steady-state oscillations 


6 « ; ; . 
Malheureusement, la démonstration donnée par CLEBSCH peut laisser place a 
quelques doutes.”’ 
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is identified as a special case of IACOVACHE’S solution. Finally, Section 5 is 
devoted to general solutions of the equations of motion in the special instances 
of rotational symmetry and of plane strain (or generalized plane stress). 

For the sake of focusing on essentials, the body forces are supposed to vanish 
throughout the succeeding analysis; the extension of the results to the case in 
which body forces are present offers no difficulties. Also, in order to avoid a 
steady repetition of tedious regularity assumptions, all such hypotheses are 
omitted from the statement of the theorems to be discussed. Suffice it to say 
that the region D under consideration need not be bounded or simply connected. 
If D is bounded, we require merely that u(x, t), together with its space and time 
derivatives of the first and second order, be continuous for (x, %), x3) in D, 
i<t<t,. If D is not bounded, u(x, ¢) must, in addition, be subjected to suitable 
regularity conditions at infinity. 


2. Duhem’s proof of Clebsch’s Theorem 


The completeness theorem announced by CLEBSCH [7] may be stated as 
follows: Let u(x,t) be a particular solution of (1) in a region of space D and for 
t<t<t,. Then there exists a scalar function p(x, t) and a vector function (x, t) 
such that u(x, t) 1s represented by (7), with 


V-p=o, (10) 


and y, \ satisfy (8) in which the wave operators 1] (n=1, 2) are given by (4). 
We now summarize, in vectorial notation, DUHEM’s [9] proof of the preceding 
theorem. Let W(x, 7?) be the Newtonian potential defined by 


ae u(é,t) dv 
Av J) Rix, é)i? (11) 
D 


W(x, t) = 


where R(x, &) is the distance between two points of D with coordinates (%1, 2, %g) 
and (&,, €,, 3), respectively. Define U(x, t) and V(x, t) by 


O= VW. Va V XW, (12) 
note from (11) that throughout D 
pei th (13) 
and recall the general identity 
VA VV.W —VxVxw. (14) 
From (12), (13), (14) follows 
a=VUPVxV, V-V=0, (15) 


which is merely the usual Helmholtz representation of the vector field w(x, ). 
Substitution for « from the first of (15) into (1), and recourse to (4), yield 


VEU elle Vv = 0. (16) 


Operating on (16) with the divergence and the curl, respectively, one obtains 
with the aid of the identity (14) (applied to V), in view of the second of (15), 


VaneU=0,0 V2? AV =0: (17) 
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JES O2U=a(x,t), O2V=b(z,4), (18) 


where 
V2a=0, V2b=0, V-b=0, (19) 


the last of (19) being a consequence of the second of (15) and the second of (18). 
Now, define A(x, ¢) and B(x, ¢) through 


tt bt 
Ant f fat a) dddt, B=, [ [ d(x, A)dadr, (20) 
a ty Oh a to 
where ¢,<fy<t,, so that 
Gel Ch aes 
et GR Lass ét? ce 
and, because of (19), 
LPA=0, VB=0,  V-B=o- (22) 
Next, let U, (x, ¢) and W(x, t) be defined by 
UsULA, bal PB. (23) 


Substituting for U, V from (23) into (18), and taking account of (21), (22), one 
arrives at 


7 Up 0 ls pee: (24) 
By (45) and (23), 
u=VU,+0xb +u*, (25) 
provided 
u* = —VA —VxB. (26) 


Equations (26), (22), together with the identity (14) applied to B, imply 

Viu*¥=0, Vxu*¥=0. (27) 
Therefore, there exists a function U,(x, ¢) such that 

u*=VU,, V2U,=0, (28) 
and the first of (28), in accordance with (25), gives 

u=VU, + VU, +-VX (29) 


Substituting from (29) into (1) and bearing in mind (4), one obtains in view of 
(24) and the second of (28) 
OPUy xs 


Ve =0, (30) 
whence 
U, = a(t) +tB (x) +(x). 61) 
By (31) and the second of (28), 
V2pB=0, V%yi=0, (32) 


Finally, define p(x, t) through 
p =U, (x, t) + U(x, t) — a(t). (33) 
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According to (33), (29), in conjunction with the second of (15), the last of (22), 
and the second of (23), 


u=VotVxb, V-b=0, (34) 
while from (24), (31), (32), and (33), 
Dig=0, LCep=o. (35) 


This completes the proof of the theorem. It should be noted that u(x, t), given 
by the first of (34), satisfies (1) provided (35) hold even if (x, ¢) is not solenoidal, 
as may be confirmed directly. 

DUHEM’s argument is easily extended to accomodate body forces and a time- 
dependent temperature field. Such an extension, in particular, assures the com- 
pleteness of a general solution to the coupled thermoelastic field equations 
employed recently by DERESIEWICcz [33]. 


3. Lamé’s Solution in the Equilibrium Case 

Since the foregoing proof of CLEBscH’s theorem does not preclude the pos- 
sibility that wis independent of the time ¢, LAME’s solution (7), (8) of the equations 
of motion (1) clearly remains complete in the special case of equilibrium. As 
pointed out in the Introduction, however, the potentials m and W are bound to 
remain time-dependent also in this event unless u(x) happens to be harmonic®. 
We now examine the manner in which LAME’s potentials involve the time when 
u is a solution of the equilibrium equations and establish the relation which 
and W bear to the PAPKOVICH-NEUBER potentials of the equilibrium theory. 

To this end we recall that PAPKOVICH’s [14] solution, which was independently 
rediscovered by NEUBER [15], has the form 


u(x) =V(O+a-W)—4(1—») B, (36) 


where « is the position vector with components (%,, %2, %3), and @(x), as well 
as W(x), are solutions of the Laplace equations 


Y27®=0, V2¥=0. (37) 


~The completeness of the equilibrium solution (36), (37) was established first in 

[14]. A more direct completeness proof was supplied later by MINDLIN [2/]®. 

We now prove the following theorem: Suppose an equilibrium solution u(x) 

of (1) is generated, in the sense of (36), (37), by the Papkovich-Neuber potentials 

D(x) and W(x). Then it is also generated, in the sense of (7), (8), by the Lamé 
potentials p(x, t) and p(x, t) defined through 


p= —2(1—-») #2 7-8 — (1-20) a-®, 
; (38) 


yp =2(1—») ue Vx¥+7xo, 


8 This requirement is merely necessary for m and to be independent of ¢ in the 
equilibrium case; conditions which are both necessary and sufficient will be given 
later. See the statement immediately following (57). 

9 See also [22] for a discussion of certain related questions. 
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provided O(x) is a particular solution of the pair of equations 
ye=4i-vy¥, V-O=2(1—»)x-®. (39) 20 


To confirm the truth of this theorem, one needs merely substitute (38) into 
(7), (8), taking account of (4), (14), (36), (37), (39), and making use of the identities 


_ 2 
= 3a—ay’ a 
V2(2-W) =a V2 +2V-¥. (41) 


With a view toward motivating the particular form of (38), (39), we observe 
that if the potentials p(x, ¢), p(x, t) and ®(x), W(x) give rise to the same solution 
u(x) of (1), they must, according to (7) and (36), satisfy the relation 


Vo t+Vp=V(P+a-P) —4(1—»)B. (42) 


Also, in view of CLEBSCH’s theorem, we may assume that (10) holds, 7.e. p may 
be restricted to be solenoidal. Next, operate on (42) with the divergence and the 
curl, respectively, and recall (8), (10), (37), as well as the general vector identities 
(14) and (41). Thus, 


ay == 419) 27-8, ey =4(1—») AVx®, (43) 


6 “Ot 
and after two successive integrations, 
hd 
Fees 416 et) : VE +t, (%) + 2(%), 
{2 
b= 2(1 =») 22 XW + ty (x) + ala). 
From (44), (10) follow 


id 


(44) 


Vb = One Vp 0, (45) 


Now operate on the first and second of (44) with (7 and (3, respectively, 
bearing in mind (4), (8), and (37). This yields 


2120) V- YE + (Vg, = V2 o, =0, 


—4(1—»Vx¥+itV?b, +7?p,=0. a 
Since (46) must hold identically with respect to ¢, we have 

V?29=0, V2o,=—2(1—2r)V-f, (47) 

Vp=0, V2b,=4(1—-») Vx. (48) 


By virtue of (41) and the second of (37), the general solution of the second of 
(47) may be written as 


Pr=—(1—2r)e-P+o, V%o,=0. (49) 


On the other hand, the second of (45) and the second of (48) together assure the 
existence of functions @(x) and A(x) such that 


~=Vx0, 7O=4(1—» ¥+7/. (50) 


10 As will be shown later, (39) always admit a simultaneous solution. 
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Substitution from (49), (50) into (44) leads to 


p=—2(1—») HOV. Y — (1-29) eB tig + oe, 


51 
= 2(1 —») Ax +h, 47x, 61) 


in which, because of (45), (47), (48), (49), and (50), 
OV "on" b,— Vw, = 0, V7O = 4(1 — 9) BV (52) 


and 4 is as yet arbitrary. 


Finally, substitute (51) into (42), use (14), (37), (52), and equate the coefficients 
of like powers of ¢ on either side of the resulting identity. In this manner one 


arrives at 
Vo, + Vx, =0, 


Vig, +V-O — A) =V[G®+2(1—v) x- WP). (53) 


Equations (52), (53) can evidently be met by setting 


A=—=4=0, m=, (54) 


if @ is a solution of (39), and (51) then become identical with (38). This accounts 
for the specific form of (38), (39). 

The pair of partial differential equations (39) for O(x) is compatible since, 
according to (41) and the second of (37), 


VP? (eV) =27-¥. (55) 


For given W(x), a solution to (39) may be constructed by observing on the basis 
of the second of (37) and the second of (39) that @(x) is biharmonic. Hence, 
in view of ALMANSI’s [23] theorem, © may be written as 


@=0'1%,0", 70'=720"=0, (56) 


in which x, may alternatively be replaced with x, x3, or x?. The two harmonic 
functions O’(x) and O”’(x), in turn, may easily be determined in such a way that 
both of (39) hold. 

The Lamé potentials (38) are seen, in general, to depend quadratically on 
the time. By (36), (37), and (41), 


V-u=—2(1—2xV-¥, Vxu=—4(1—»Vx. (57) 


Consequently, p and p in (38) ave independent of the time in the.equilibrium case 
if and only if the displacement field is both equivolununal and irrotational, which, 
by virtue of (1), requires that u(x) be harmonic. As is apparent from the theorem 
at the beginning of this section, every equilibrium solution u(x) of (1) may be 
represented as a sum of an irrotational and an equivoluminal motion the time- 
dependence of which cancels upon addition. This decomposition is evidently not 
unique. 

To illustrate the foregoing conclusions we turn to a specific example. Thus, 
consider KELVIN’s elastostatic solution! appropriate to a concentrated load at 


11 See LOVE [25],. art. 130. 
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a point of a medium occupying the entire space. Let the load have unit magnitude 
and be applied at the origin in the direction of the x,-axis (a= 1,2, 3\ug dane 
corresponding PAPKOVICH-NEUBER potentials then take the indicial form 


Oxi 
P=0, t= 162(1—v) wR’ (58) 


where R is the distance from the origin, while 6,; denotes the Kronecker delta. 
In accordance with (36), the equilibrium displacement field generated by the 
harmonic stress functions (58) has the components 


a fo ats a als 
Ui; 162(1—v) wR R2 a (3 49) at (59) 

In this instance one finds that (39) may be met by taking 

go (60) 


so that (38) yield the associated Lamé potentials 


en 2(1—») we? {+2 
= aa — Le. 
emi anak | oR? | ; (61) 
YW. = Saif Hj | we = 
i eaiiae || lar 


in which e,,;; are the components of the usual alternator and summation over 
repeated indices is implied. The stress functions (61) conform to the wave 
equations (8). Using (61) in conjunction with (7), we arrive at the decomposition 


u, (x) =u, (x,t) + uj (x, t), (62) 
where 
p_ -OGj fl 2(1—v) wt? 3% ey %j De ae Hy Xi 
men inte oR? ( R bai) Sete, ee éai)]) (63) 
Ae ep Ow; a 1 \ Hee (0 = 3a) ay eee Key Xj 3 
4 tkq OX, ~ Sauk o R? at R2 T x4 R2 . 


The irrotational and the equivoluminal motions represented by w;(x,?) and 
u; (x,t), respectively, individually satisfy the equations of motion (1). 

The resolution of an equilibrium solution into dynamic component solutions 
has an analogue in the theory of the biharmonic equation and of the repeated 
wave equation. To make this clear we recall a result established separately in 
[19] and contained in a more general theorem due to Bocaio [24]: every (suf- 
ficiently regular) solution g(x, ¢) of the equation 


LiiCeg=0 (64) 
admits the representation 

9(%, 1) = G9 hg, (65) 

where the functions g’(x,t) and g’’(x, 2) satisfy the ordinary wave equations 

Oig’=0, Og’ =o. (66) 


The operators (];, [12 appearing here are again given by the first of (4), but 
c; and cy, at present may have any distinct non-zero values. 
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If, in particular, g happens to be independent of ¢, (64) degenerates into the 
biharmonic equation V4 g=0, the complete integral of which admits ALMANSI’S 
[23] representation in terms of two harmonic functions. Thus, 


g(x) =h'+2-h’", (67) 
in which h’(x) and h’’(x) satisfy 
V?eh’=0, Veh" =0. (68) 


The representations (65), (66) and (67), (68) may be regarded as respective 
counterparts of LAME’s general solution (7), (8) to the equations of motion and 
of the Papkovich-Neuber solution (36), (37) to the equilibrium equations. Boccio’s 
theorem, and hence (65), (66), still hold if g is a function of position alone, but 
in this event the resolving functions g’ and g” ordinarily remain time-dependent. 
For instance, the scalar biharmonic function 


g(x) = x4 (69) 
is the sum of the two wave functions 


»__ (ai +c7t) n+) 
& n> ck—c2 ra & ce—c3 (70) 
4. Connection between Lamé’s and Iacovache’s Solution 
IACOVACHE’S [16] solution of the equations of motion (1) is given by 
u(x,t) =2(1—7) 07 G—VV.-G, (71) 
where 
fel? 5G 0: (72) 
Our present aim is to examine the relation between (71), (72) and LAME’s solution 
(34), (35). 


Invoking once again the special case of Boccio’s [24] theorem arrived at 
independently in [19], we note first that there exist functions G’(x, ft), G’’(x, t) 
such that 


G(4,ij)=G 4G’ (73) 
and 
HiG’=0, OG’ =o. (74) 
From (4) and the second of (74) follows 
21-7) DIG’ =77a@". (75) 
Using (73), (75), we may write (71) as 
u=V? GY — VV -(G! + G’) (76) 


from which, on applying the general identity (14) to the vector field G’’(x, 2), 
ciau, u=—VV-G! —VxVx@". (77) 
If, finally, we define functions p(x, ¢) and p(x, ¢) by 

g=—V-G, b=—VXxe", (78) 
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(71) takes the form 

u(x,t) =Veo+tVxp, V-p=o, (79) 
while (74), (78) yield 
Dig=0, DCep=o. (80) 


Equations (79), (80), however, are identical with (34), (35). LACOVACHE’S solution 
thus reduces to Lamk’s solution of (1) if the stress function G(x, t) 1s subjected to 
the transformation (73), (74), (78). In view of the completeness of [ACOVACHE’S 
solution (71), (72), established previously in [19], the preceding reduction supplies 
an alternative proof for the completeness of LAmg’s solution 34). (35) 5 im 
comparing this proof with DuHEm’s [9] argument, given in Section 2, it should 
be remembered that we have at present assumed the generality of IACOVACHE’S 
solution which was proved in [79] with the aid of the theory of retarded potentials. 


If G is independent of the time, (71), (72) degenerate into GALERKIN’S [18] 
general equilibrium solution 


u(x) =21 NV G—RV-G, (81) 
V*G=o, (82) 


of which IACOVACHE’s solution is therefore a dynamic generalization. GALERKIN’S 
solution (81), (82), in turn, as pointed out by MINDLIN [2/], is reducible to the 
PAPKOVICH-NEUBER solution (36), (37) by setting 


Oa)=—V.G—a2-¥, Wj 177 e: (83) 


Similarly, as observed in [19], the time-dependent transformation 


O(x,)=—V-G=e2% Waj)=]—207G, (84) 

applied to G(x, t) of (71), (72), leads to 
u(x,t) =V(®+a-¥P)—4(1—7r) BP, (85) 
Di®=]=—2: fi, bade eee (86) 


which may be regarded as a dynamic generalization of (36), (37). Since the 
potentials ® and Y in (86) are coupled, this form of the complete solution to 
the equations of motion (4) is of no practical interest. 

In attempting to judge the comparative merits of [ACOVACHE’s solution (74), 
(72) and of Lame&’s solution (79), (80), the following considerations would appear 
to be pertinent. IAcovACHE’s solution remains complete in the equilibrium case 
if the generating potential G(x, t) is taken to be a function of position alone. 
In contrast, as explained earlier, LAm#’s solution has no equilibrium counterpart 
in the foregoing sense. This circumstance reflects the relative economy of LAME’s 
solution to the equations of motion, which is simpler in structure than the available 
complete solutions to the equations of equilibrium. It is clear from (73), (74) 
that the single stress function Gin (71), which satisfies the repeated wave equation 
(72), is equivalent to two stress functions each of which meets an ordinary wave 
equation. The number of independent solutions of the wave equation involved 
in (71), (72) and (79), (80), respectively, is thus the same. On the other hand, 


Potentials for Elastic Vibrations 45 


(71) contains space derivatives of G up to the second order while only first- 
order derivatives of y and w are seen to enter (79). Finally, (79), (80) are con- 
veniently transformed into general orthogonal curvilinear coordinates, whereas 
(71), (72) give rise to exceedingly cumbersome forms when referred to curvi- 
linear coordinates, with the exception of circular cylindrical coordinates. For 
all of these reasons, LaM&’s solution deserves preference over IACOVACHE’s general 
solution to the equations of motion in applications of elastokinetic theory. 
NoLt [20] recently dealt with the integration of the equations of motion in 
the case of steady-state oscillations of elastic bodies. He thus considered motions 
of the type 
THOM AIS Ca Weed (87)? 


and showed that if u(x, 7?) given by (87) satisfies (1), the amplitude function 
Uy (x) admits the representation 


Uy (x) =2(1 —») V2, + (1 — 2”) ~ Gee.G. (88) 
w? w? 
Vat =| if ! | Go= 0. (89) 


The solution (87), (88), (89) is identified as a special case of IACOVACHE’S solution 
by taking 

Got) =Golx) e°?" (90) 
in (71), (72). A more useful complete representation of the steady-state vibration 
(87), however, is obtained by setting 


p(%,t) = po(x)e?, (x,t) =ho(x) e'” (91) 


in LAME’s general solution (79), (80). 


5. Axisymmetric and Two-dimensional Solutions of the Equations of Motion 
Let (7, y, z) be circular cylindrical coordinates defined by the mapping 
My=7COSY, %,=7siMy, %%=2, 


(92) 
Oise Core Oa = DT) IO 2 <0, 


and let (u,,,,w,) stand for the corresponding cylindrical components of the 
displacement vector. A displacement field possesses (torsionless) rotational sym- 
metry about the z-axis if 


U,=U,(7,2,t), U,=U,(7,2,t), m=O. (93) 


Suppose we denote the cylindrical components of the vector potential in 
Lamg&’s solution by (p,, p,, y,). Evidently, (79), (80) represent a solution of (14) 
of the form (93) provided we choose 


9=9(7,2,1), Ww=¥24), YY =¥.=09, (94) 


12 The fact that w here designates the circular frequency while w in (5) denotes 
the rotation vector ought not to cause confusion. 
13 2(1— v) Gy=q(z) in the notation of [20]. 
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and the associated non-vanishing cylindrical displacement components then 


become 


og ee Op i Oey 95 
La Tye Cnpig te Ms az | ape Fy (95) 


The two scalar potentials p, p satisfy the wave equations 


ig=0, Leyp=o, (96) 
in which (]2, (33 retain their previous definition (4), and the Laplacian operator 
now appears as ge on 22 : 

VS Gar ae (97) 


Equations (95), (96) constitute the axisymmetric form of LAMEe’s solution. 
A rotationally symmetric version of IACOVACHE’s solution is obtained from 


(71), (72) by setting G=GH,2h, CSE (98) 


where (G,, G2, G3) are the cartesian components of G. Consequently, 


Ip Ut, = 2A = 9) VG — 


if Orv 02? ck 


Oz? ( 
ERE =0; (100) 


and (1%, (12 have the same meaning as in (96). If, in particular, G here is inde- 
pendent of the time, the foregoing solution coincides (except for a constant factor) 
with Love’s [2 ]!4 general solution of the equations of equilibrium in the presence 
of rotational symmetry. 

The connection between the alternative axisymmetric solutions (95), (96) and 
(99), (100) follows at once from an appropriate specialization of the general 
conclusions reached in the preceding section of this paper. Thus, (73), (74), (78), 
because of (94), (98), yield in the present instance 


G(r, z,t) =G'+G", (101) 

DiG’=0, [G6"=0, (102) 
CoG’ OG 

CS a Co) 


The transformation (101), (102), (103) reduces (99), (100) to (95), (96), as was 
observed independently by PREDELEANU [26]. Equations (95), (96) and (99), 
(100) individually represent complete solutions of the equations of motion (1) 
in the rotationally symmetric case governed by (93). This is readily shown to 
be true of IAcovACHE’s solution (99), (100) on the basis of a suitable adaptation 
of the general completeness proof presented in [19]; the completeness of Lamt’s 
solution (95), (96) may then be inferred directly with the aid of (101), (102), (103). 

We tum, finally, to a discussion of general solutions of the equations of motion 
in the two-dimensional theory of elasticity. In the event of plane strain (parallel 
to the plane x,=0) the displacement field is characterized by 


Uy = 4 (%y, %y,t), Wy = Ug (%1, Xp, t), ug =O. (104) 


Bewee att. 183) po 274 
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It is clear from the two-dimensional counterpart of CLEBSCH’s theorem (Section 2) 
that Lameé’s solution remains complete in the present instance if the generating 
potentials are taken in the restricted form 


P=P(%,%,t), Ye=Y(%,%2,t), Yr=Y.=O0. (105) 
Equations (79), (80) now lead to 
er as OP BEE ge 
ira nir oye caart oz: to) 
197=0, Cey=o. (107) 


A plane-strain version of IACOVACHE’s solution, on the other hand, is obtained 
by setting 
in (71), (72). This choice, after some elementary manipulations involving (4), 
results in 


a?G 2G (1—2r)0 @G 


Co? O2G=0. (110) 


Equations (109), (410) are identical (except for a constant factor) with the 
complete two-dimensional solution deduced in a different manner by SOBRERO 
[27]. SOBRERO’S general solution is reducible to the two-dimensional form (106), 
(107) of LAmg’s solution by employing once again the transformation (73), (74), 
(78). Bearing in mind (105), (108), we thus arrive at 


Gi eG Gr (414) 
ei, © eel OP Re (112) 

mS Peace her hoe” 
? ee ee (113) 


which link the two solutions (106), (107) and (109), (110). A plane-strain solution 
of (1) in terms of two scalar potentials each satisfying the iterated wave equation 
(110), was given by G. C. Motsir [28]®. Motsi’s result is also found to be a 
special case of IACOVACHE’S solution (71), (72) and is reached (except for a constant 
factor) by choosing 


Gy = Gy (%1, %2,t), Gy=G_(%,%,t), G,=0. (114) 


The present discussion has so far been confined to plane strain. As is well 
known, the transition to the corresponding results for generalized plane stress 
is effected by replacing Porsson’s ratio v with »/(1-++-v). While the wave velocity 
C) in (4) remains invariant under this change of Polsson’s ratio, ci passes over 
into ¢?, where 


NOS 2 
2 Eee: ne 


15 See Equations (116) on p. 384 of [28]. The coefficient 442m appearing in 
both of these equations is in error and should read A+ yp. 
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Equations (106), (107) have been employed in a number of recent two- 
dimensional studies dealing with uniformly moving loads, punches, cracks, and 
dislocations! as, for example, by SNEDDON [30]. In contrast, RaDoK [31] 
reconsidered several plane problems of motion on the basis of a dynamic analogue 
of Arry’s elastostatic stress function. The same scheme was encountered later 
on by TEODORESCU [32]. 

The comments made at the end of Section 4 concerning the advantage of 
Lamk’s over IACOVACHE’S form of the general solution to the three-dimensional 
equations of motion, appear to be equally valid in connection with axisymmetric 
and plane problems. In fact, it is probably safe to say that the various displace- 
ment potentials (or stress functions) which have been introduced into elasto- 
kinetics as a result of efforts to generalize known stress functions of elasto- 
statics, while interesting from the theoretical viewpoint, are of limited interest 
from the point of view of applications. For the requirement that a complete 
solution of the equations of motion reduce to a complete solution of the equi- 
librium equations whenever the generating stress functions are independent of 
the time precludes the utilization of an element of relative simplicity which is 
inherent in the dynamic case but absent from the degenerate equilibrium case. 
This fact is reflected clearly in the comparative economy of LAME’s classical 
general integral of the equations of motion which has no such counterpart in the 
equilibrium theory. 


Note Added in Proof. The author recently encountered certain publications of 
which he was previously unaware and which necessitate the following additional 
comments: 


(1) Correct proofs of CLEBSCH’s [7] completeness theorem were also given inde- 
pendently by Somici1ana [34]!’ (1892) and by TEpoNeE [35] (1897) prior to the 
appearance of DUHEM’s [9] proof in 1898. Also, it is evident from the introductory 
remarks in [34] that SomiGLrana recognized the shortcomings of CLEBSCH’s earlier 
argument}!’. Nevertheless, DuUHEM’s proof, which is summarized in Section 2 of the 
present paper, appears to be more economical and direct than the alternative proofs 
contained in [34], [35]. 


(2) Further, SomicLiana [34] (1892), deduces the representation!” (71), (72), from 
Lame’s solution (7), (8) in a manner which assures the completeness of (71), (72). 
Consequently, the general solution (71), (72), which here and throughout the con- 
temporary literature has been attributed to IacovacHE [16] (1949), should instead 
be credited to SOMIGLIANA. 


(3) Finally, reference should be made to recent papers by ARZHANIKH [36] and 
by Sotomon [37], which contain variants of the solution (85), (86) in terms of stress 


functions satisfying a pair of interlocked partial differential equations of the second 
order. 
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Korn Inequalities for the Sphere and Circle 


B. BERNSTEIN & R. A. TOUPIN 


1. Introduction * 


Let © be a given region in either two or three dimensions. Let ¢” denote 
the class of all vector fields u(a), continuous in the closure G of & and 2 times 
continuously differentiable in © such that the Dirichlet integral 


6 


is finite. Let d;; and w,; denote the symmetric and antisymmetric parts of the 
gradient of wu: 


di; = Uy 57) = 3 (4%, + %,5), (1.2) 
C055 = Mes, 5) =F (M;,49 — 4,2) - 
We shall also use the abbreviated notation 
a —d;,d;;, wo? = 0;;0;;; Pid a re (1.3) 


The summation convention, Cartesian coordinates, and Cartesian tensor 
notation will be used throughout. 

A Korn inequality states the existence of a number K which depends only 
on the shape of the region & such that 


fordvsk f ddv (1.4) 
6 6 


for all vector fields which satisfy one or more side conditions. That an inequality 
of the form (1.4) cannot hold for arbitrary vector fields in @© follows from the 
existence of pure rotations, u,=¢;;, %; ,, h =constant vector, ¢93=1, 6; ,= i; A] 


* While this paper was being written, the authors were in touch with L. E. Payne 
& H.F. WEINBERGER, who, by using a variational approach, were able to obtain 
extrema for the expression (1.6) in the second case for the sphere and the circle. 
But they were not at that time able to prove that their largest extremum was actually 
the supremum of (1.6) and hence Korn’s constant in the second case. However, 
after this paper had gone to press, the authors were delighted to learn that PayNnE & 
WEINBERGER had succeeded in formulating a proof of this, and thus Korn’s constant 
in the second case is now known to be 43/13 for the sphere and 3 for the circle. 
Payne & WEINBERGER will publish their work in a subsequent issue of this Archive. 

The authors were dissuaded from withdrawing this paper in the light of the 
results of PAYNE & WEINBERGER because, although the numerical values of upper 
bounds for Korn’s constant contained herein are no longer of value, the remainder 
of the work remains a contribution. Also the adaptation of the method of FRIEDRICHS 
contained herein is applicable to more general regions, whereas PAYNE & WEINBERGER 
presently feel that their method is not. 


4* 
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We define, with Frrepricus [1], Z® to be that subset of @* consisting of 
those fields ‘w which vanish outside some compact subset of &, and Z° to be 
the set of fields win @! for which fields *u in Z® may be found to make D(u—‘u) 
arbitrarily small. 

Korn inequalities for a class of regions, called {2-domains by FRIEDRICHs [J], 
are known to exist. For a definition of the class of 2-domains the reader is re- 
ferred to [Z]. It suffices here merely to note that this class of regions includes 
the sphere and the circle. Moreover, for an Q-domain, Korn inequalities exist 
under three types of side conditions. Following Korn and FRIEDRICHS, we 
consider separately the three sets of side conditions which, using the terminology 
of FRIEDRICHS, we call the three cases: 


First case: UucD®, 

Second case: fw,dv=0, uc, (1.5) 
6 

Main Case: fa ;dv=0, uc, and d;,;=0. 


If a Korn inequality exists for a given region &, subject to a given one of 
the side conditions (1.5), then the set of real numbers which are the values of 


fa dv 
10) 


cae (1.6) 
® 


subject to the given side conditions, has an upper bound, and hence a least 
upper bound. We call this least upper bound Korn’s constant for the region & 
an the first case, the second case, and the main case respectively. We denote these 
dimensionless positive numbers by K&, KG, and KX. Any number greater than 
or equal to KG, Q=1, 2, M, will be called an upper bound for K@, and any 
number smaller than or equal to A§ will be called a lower bound for K&. Sub- 
stitution of any upper bound for a Korn constant in (1.4) yields an inequality 
valid for all vector fields satisfying the appropriate set of side conditions (1.5). 

The proof of the existence of a Korn inequality in the first case for an Q- 
domain follows immediately from the identity 


JL Oh died de dor alma le das (37) 


which is shown by FRIEDRICHS [1] to hold for all win D° when © is an Q-domain. 
Clearly, then, for the first case 


SRB SINS (4.8) 


and Kos! for all ©. Moreover, since it is easy to construct vector fields Uc @g® 
for which trd=0, we see from (1.7) that K6=1 for any Q-domain. 


‘The useful role of the identity (1.7) and of Korn’s inequality in the first 
case which follows from it in studies of the displacement boundary value problem 


of linear elasticity theory is made apparent in the papers of Korn [2, 3] and 
FRIEDRICHS [1]. 
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To prove the existence of a Korn inequality in the second and main cases 
is not so easy. As shown by FRIEDRICHS and in Section 4, the second case follows 
easily from the other two, whence the name, the “‘main case’. We base our 
understanding of the main case on the paper by FRIEDRICHS. With him we 
confess inability to follow Korn’s original treatment. 

The constructive proof in [1] of the existence of a Korn inequality in the 
second and main cases provides a technique for computing, at least in principle, 
upper bounds for the Korn constants K¥ and K% for any Q-domain. But to 
our knowledge, no one has found a value for the Korn constants K% or KY 
for any region. Nor have we seen in the literature any upper bound for these 
numbers for any region whatsoever. Our curiosity was aroused as to just how 
large or how small the Korn constants might be for even a region so simple 
as a sphere or a circle. For these regions it turns out that FRIEDRICHS’S construc- 
tion for the second and main cases can be simplified somewhat, though we retain 
its essential ingredients. To the simplified proof for the sphere and circle we 
add a determination of the following upper and lower bounds: 


For the sphere: 3S .K4<20, 35.K¥ <17 


A 
Por thecucle: "2 = Keo eA! (1.9) 


It is our belief that the method developed below for the sphere and the 
circle can be extended to more general regions. 


2. Lower bounds for Korn’s constants 
a) The case of three dimensions. For the purposes of this section let © be 
any region for which the centroid and moment of inertia tensor about the centroid 
(EULER’S tensor) have meaning. Let the origin of the Cartesian coordinate system 
x, coincide with the centroid of &, so that the moment of inertia tensor about 
the centroid is given by 


Lim [ox 8,405 (2.1) 

6 
Let h be a constant unit vector field, 4,4; =1, and consider the vector field 
Ones (7, Xn) Ciba hy Xq: (272) 


Viewed as a displacement field of an elastic body say, the field wu represents 
a uniform twisting in which each particle in a plane normal to h remains in 
that plane but is displaced in a direction normal to h and its position vector «x. 


We obtain for the vector field (2.2) 
eee == = (0;; Fz h, h;) I;;, 
fwrdv —— fd?dv + 2h; h, I; 


72 


6 6 
q;; j = O ) 
from which it follows that 
2d 
ee, ATs hic hy 


aed | oie) Tey (2.3) 
G 
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Also, since the derivatives of u, are linear and homogeneous in x,, and since 
the origin is at the centroid of &, (1.5)2 is satisfied. Furthermore, since d;; ;=0, 
our example belongs to the main case. 

In general, the right-hand side of (2.3) will vary with the choice of the direc- 
tion h. As can be shown by the method of Lagrange multipliers, it will have 
extrema when h is in the direction of the principal axes of the moment of inertia 
tensor for the region. The largest we can make the right-hand side by varying 


the unit vector h is id eee (2.4) 


where Imax =Imea=Imin are the eigenvalues of the moment of inertia tensor. 

Considering different regions now, (2.4) will have the smallest value for a 
region for which the three eigenvalues are all equal, in which case its value 
will be 3. Thus for any region in three dimensions K@2>3, ke 23. Smee (24 
yields a lower bound for K¢ (and KG) for any given region we may make the 
following observations: 

4. There is no single number K which will yield a Korn inequality in the 
second or main cases for all three-dimensional regions—z.e., the set of Korn 
constants in the second and main cases for three-dimensional regions has no 
upper bound. 

2. Given any region © in three dimensions, one can find a sufficiently long 
and thin right circular cylinder &’ whose Korn constant KY, exceeds KY. 

b) The case of two dimensions. Next we establish a lower bound for Korn’s 
constant K§ for regions in the plane. For this purpose we compute the ratio 
(1.7) for the 2-dimensional vector field w having the form 


U; = (hy, Xp) €55 %;- (2:5) 


All indices now take the values 1 and 2, and e,,=—1, euij;)=%,- Again we assume 
the origin of the coordinate system coincides with the centroid of the planar 


region and let iy =I x, dv (2.6) 


denote the moment of inertia tensor of the planar region ® about its centroid. 
For a vector field having the form (2.5) we get 
2 dit — Ji, h;h;, 
6 
2fordv=4],h,h,, 
6 
from which it follows that 
f adv 
© 4 Jig hy hj 
[Pav ~ rt J— Tyhihy | (2.7) 


The maximum value of the right-hand side of (2.7) is attained when h is an 
eigenvector of j;; corresponding to its maximum eigenvalue. Thus the maximum 
value for the ratio in (2.7) is Ano 

Jmax+2 Jmin , (2.8) 
which lies between # and 4. It follows that for any region in the plane, K2>4 
The lower bounds announced in Section 1 are established. “ae 
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3. Some differential identities and algebraic inequalities 


In what follows we shall make repeated use of the identities 


Ojn,1= 4; — TnI; (3.4) 
Gig nat Gnt,33 — Un,g1— rex = 05 (3.2) 

and of the inequality, 
(trd)?< nd?, (3.3) 


where m is the dimension of the space. The differential identities (3.1) and (3.2) 
are well known in elasticity theory. The latter are the familiar compatibility 
conditions that d;; be expressible in the form (4.2),. Inequality (3.3) follows 
easily from CAUCHY’s inequality, 2ab<a?+}?, in a frame of reference where 
d;,is diagonal. Since both members of the inequality are invariant under ortho- 
gonal transformations, (3.3) is valid always. 


From (3.1) we get 
O55, ij, = 2 (Giz Gig — Gig Tei,j) - (3.4) 
We shall be concerned in particular with vector fields which satisfy the 
differential equations d;; ;=0 of the main case. For this class of vector fields 
we get from (3.1) and (3.2): 
V?0;;=@;;,45=0, 
Oi5,5 = — 4ij,a) (3.5) 


G5 uk a ds k,i7- 


4. Upper bounds for K% in terms of an upper bound for Ke 


That every upper bound for Korn’s constant in the main case determines 
a corresponding upper bound for Korn’s constant in the second case was shown 
for an Q-domain by FRIEDRICHS [1]. We wish to present here a slight variation 
or refinement of his result which is necessary for us to establish upper bounds 
for K%, as small as those announced in Section 14. 

In this section we take & to be an 2-domain to which the divergence theorem 
applies. As shown by FRIEDRICHS, any vector field in @1, satisfying the side 
conditions (1.5), of the second case, may be uniquely decomposed into a sum 


“a= Ut”u, (4.1) 
where ‘u€QG°, and "wu satisfies the side conditions of the main case, namely, 
"Uc3, TT Sp 0, f'@,; dv ==), (4.2) 

6 


where 'd;; and 'w;; are the deformation and rotation measures of “u. This 
follows from the known uniqueness, existence, and differentiability properties 
of solutions of (4.2) with prescribed data for ‘’“w on the boundary, and the fact 
that f{ ’w;,;dv=0 for all ’u. 

6 


Let K” be an upper bound for K%, so that we have 


Sorde =< K" f "@dv. (4.3) 
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Now given a field w€@1, decomposed as in (4.4), let *u be any field in 7 Z® 


and define 
*U — "EY +- Moy) . (4.4) 


Then 
f*w2dv aa +o? + 2°0,; "0, ,] dv 
3) 


=f [°a? a (tr *d)? aE je "1 2 + 2°U; ; "0; ; | dv. 
10) 


(4.5) 


In (4.5) we have used (1.7) to eliminate ‘w® and (4.3) to eliminate “w*. Using 
the identity (3.5), and the fact that ‘w vanishes outside some compact subset 
of ®, we can transform the last term in (4.5), with the aid of the divergence 
theorem, as follows: 


2 f “Uj "Oj 5 dv = Dp if [(°%4; "O; ie == “U; "Oi55] dv 
6 1) 


= 2h td dv = — 2f (ted) (ted) dv (4.6) 
103) 


Ce HRC 
< f [(tr°d)? + (tr “d)?] dv, 
6 
where the last line follows from CAUcHY’s inequality. Substituting (4.6) into 
(4.5) and applying (3.3) yields 
ftwrdvs fiddv+(K” +n) f["'ddv. (4.7) 
6 6 6 
FRIEDRICHS [J] has shown that 
f *d?dv= f [d? +d" dv, (4.8) 
6 6 


and thus (4.7) and (4.8) yield 

f *o?@dvu<(K" +n) f *d2dv. (4.9) 

6 6 

Now, as remarked by FRIEDRICHS, any “uC G® as well as its first derivatives 
may be approximated uniformly by a *u€ Z*M¢?. Thus (4.9) holds for all fields 
(4.4), where *u€Z*. Moreover, since ‘w€D°, we may find a ‘we€ Y*® to make 
D(*u—wu)=D(*u—'u) as small as desired, and since 
D(u) =f @dv+fordv 
6 Gs 

for any u, it follows from (4.9) that 


fordv<(K"” +n) fa@dv (4.10) 
6 6 
for any u satisfying (1.5).. 
Therefore if K’’ is an upper bound for Korn’s constant in the main case, 


then K given by eae, (4.11) 


is an upper bound for Korn’s constant in the second case. 


FRIEDRICHS showed that if K’ is an upper bound in the main case, then 
K given b 
6 y K = 2max(1, K”) (4.12) 


is an upper bound in the second case. 
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Since we have established that, in three dimensions, K’”’=>3, the upper bound 
for the second case as given by (4.11) is never worse than the upper bound (4.12) 
for any allowed value of K” if the dimension of the space be three. In two 
dimensions, since we do not know if K’’>2, we are unable to make a similar 
comparison. 


5. The main case for the sphere and the circle 


In this section we consider only vector fields u satisfying the differential 
equation d;; ;=0 of the main case. 

Our proof of a Korn inequality in the main case and our calculation of 
upper bounds for the Korn constants K% for the sphere and the circle rest ulti- 
mately on the construction of a function y defined in the sphere G (or circle ©) 
of radius R with the following properties: 

1° y>O in © (in ©). 

2° p and its first derivatives are continuous in & (in ©). 

3° There exists a sphere ©, (a circle ©) concentric with S (with ©) of radius 
Rs<R such that 

a) In ©; (in 5) y is equal to a function which is twice continuously differen- 
tiable in ©, (in ©). 

b) In G@—G, (in © —@;) p is equal to a function which is twice continuously 
differentiable in the closure of G— ©; (of €—&;). 

4° V2 w=h(x) in © (in ©) where h=1 in G@—G, (in ©—@,) and h=— ¢, 
in ©; (in @;) where c, is a positive constant. 

5° p and its first derivatives vanish on G, the boundary of © (on C, the 
boundary of ©). 

By GREEN’s theorem and properties 4° and 5°, it will follow that the mean 
value of h(a) over © (over &) must be zero; hence, the constant c, is not at 
our disposal but is fixed by the radius R,. 


For the sphere yw is given by 


yw=i[7+2R3/7—3R*] in 6—G,, (5.4) 

yp =o l(R3 — RB) P/R5 +3. R?(R — R,)/Re] in S. (5.2) 
For the circle y is given by 

y=1[7?— R?+ R*1n(R4/r*)] in ©—@,, (5.3) 

yp = 2 [(R5 — R?) P/R5 + R? ln (R*/R)]_ in Gs. (5.4) 


For the remainder of this section, and in Section 6 until (6.18) as well as in 
Section 7 until (7.9), we add to the requirements (1.5)3 the restriction that the 
first and second derivatives of u; be continuous in © (in ©). We remove this 
restriction after (6.18) and after (7.9). 

The objective of the remaining calculations of this section is to obtain in- 
equality (5.19), which we shall use in the following two sections to obtain Korn 
inequalities for the sphere and circle. We shall understand below that, although 
we ostensibly are treating the sphere, we may, in the rest of the calculations 


of this section, replace ©, G,, etc., by ©, G5, etc. 
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Although any choice of R, between zero and k will yield functions yw for 
(5.4) and (5.2) or for (5.3) and (5.4), we shall leave the choice of R; arbitrary 
in this section, for not until Sections 6 and 7 will it be apparent how R; may 
be most advantageously chosen. 

Although we are dealing with the main case, not until after equation (5.15) 
shall we assume that the mean value of the vorticity (w,;) vanishes, and until 
then our relations will be independent of this condition, though they will depend 
on the condition d;; ; = 0. 

From 4° we have 


f wdv=cfo%dv+ fhwrdv. (5.5) 
S—S5 S6 S 


By repeated use of the divergence theorem, dropping boundary terms because 
of 5°, and by using (3.5),, we get 


ees =a Pr nr Oj Dig TY = 2 J Os, 4 O52. (5.6) 


Then adding f w?dv to both sides of (5.5) and using (5.6), we obtain 
Ss 


Jo®dv= (1+) for'dv+2 f yo,;,0;; ,dv. (5.7) 
S Se S : é 
Using PoINCARE’s inequality [4] in the form 


fv SES 05,055,040 +o" Va, (5.8) 


where P;>0 depends only on the region G;, V; is the volume of ©,, and 


we get from (5.7) 
Busey ei ae 
J v= (1%) Py J Oi, Oi, ne ac 2 J Pej, ni, dv+(1t4a4)V O°. (5.10) 


Now since y>0 in iti ini in S 
y> S, p has a positive minimum in &,, say c,R?. Then 


1< 5g in on 
and hence 
2 
i dv S[(1+ 4) Pilce Ro] J vei Ory dv + a PO; n Oy 2 dV + 


cacy), @ Soy J por; ,0,;,.40 + (14+) Vso, ay 
where : ; 
C= (1+) Bie, R3+2. (5.12) 
From (5.11) and (3.4) follows 


2 <a 7 as 
J? dvs 205 Jp (dis nai, — 45,4 Ey5,5) dv + @” V5(1+ cy). (5.13) 


Korn inequalities 59 


Using the divergence theorem, remembering that y=O on GS, and d;; ;=0, 
we get from (5.13) 


Ba 725, | _ SV, 0 (4?) 2 — dig nn dig +¥, 2 Gj Mei | dv + @? Vo (1+ G1). (5.14) 


Next we use the identity (3.5), and the fact that both y and y , vanish on 6 
to obtain from (5.14) . 


Lee & of [Vand + 29 67 (bun Gig — 4_4y;)] 0 + @*Ve(1 +c). (5.15) 


Now suppose that w also satisfies the second side condition (1.5), of the 
main case. If we define 


U; ps x; Digs 
we shall have 
"hg ip see rag Ori, (5.16) 
i *0,;, dv = O : 


Thus *d;; ;=0, and we may substitute (5.16) into (5.15) to get 
Zh *or dus oJ [pend + 2p 5; (dyn dj; — 4jp4,;)| dv. (5.17) 
But since (1.5)3 holds for @,,;, 
REE OA SITN A ASC (5.18) 
where V is the volume of 6. 
Thus from (5.17) and (5.18) we conclude that 


f° dvs oJ [wand + 29 15 (den dig — Gin dn,)] Er, (5.19) 


which is the relation we sought in this section. An inequality of the form (5.19) 
holds also for the circle. 


6. An upper bound for Korn’s constant for the sphere in the main case 


We proceed now to use the results of the previous sections to obtain an upper 
bound for Korn’s constant for the sphere in the main case. 


To simplify the notation we set 


A= Rik, 7 Sale (6.1) 
For the function w defined in (5.1) and (5.2) we have 
Y 7 = 5 10,; 1 — RAP) + 3 (Be?) nm, 2,) in S— Ss (6.2) 
where 
N= x7, nn,=1, (6.3) 
and 
Ae 1) 0;; In Ss (6 4) 
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Now the minimum value of y in G;, which is, by definition, c,R3, is assumed 


on &;, and is given by 


oa 
gma eR. (6.6) 


Thus c3 as given in (5.12) has the value 


6/3 P; 
C= 37-3041 RE 7 io) 


Insertion of (6.3), (6.4), and (6.5) into (5.19) yields 


fordvs f[—q@—a2(trd)") dv+2 f{ |d24 2(trd)24+ 2 A Aijnny| d v, (6.8) 
S Ss S-S;5 
where 

A,, = — (tr)? 6,; +d? 6; + 3 den di; — 347 Te:- (6.9) 
Since the first integral on the right-hand side of (6.8) is negative, we may drop 
it to obtain, using inequality (3.3) in the second integrand, 


fordv<? f [72+ 724i; npn]] do. (6.10) 
S S—GS6 


Consider now the form 4; ;n,; ,; at one point in S — ©;. If there we diagonalize 


d;;, getting diagonal Soren ae d,, d,, d,, we have, remembering that n;n,=1, 


4 


Aj; 1; 2; = d, do(1 — 3g) + dy d3(1 — 3mt) + dy d,(1 — 3.03). (6.11) 


IEA 


By the method of Lagrange multipliers, we see that the extrema of the 
form (6.11) are attained when n,; has the form (1,0, 0), (0,14, 0), or (0, 0, 4). 
In the first case, using CAUCHY’Ss inequality, we get 


A; N,N; = — 2d, d, + dy dz + dy d,< dj + d3? + 
_GtG , Brat _ 3di+3di+2d} _ 3 a (6.12) 
] 2 | 2 5 S15 : 


Inequalities similar to (6.12) may be obtained for the other extrema of (6.11), 
so that for all 1; we have 
A;; n,n; S 2 d?. (6.13) 


4p Eo —— 


But (6.13) is invariant to orthogonal transformations of the coordinates, and 
thus it holds for all Cartesian systems. 


Putting (6.13) into (6.10) yields 


Let 2 f 7+3 adv. 
* fr+s%) or 
Now in S — ©, we have R3)73 <3, so that from (6.14) we may conclude that 
C. 
fordos? (743%) [dvs (74329 farar. (6.15) 
8 S 


S S-S 
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A Korn inequality for the sphere in the main case is proven under the restriction 
that the first and second derivatives of uw; are continuous in G. There remains 
the problem of determining the best, 7.e., the smallest value of the coefficient 
in (6.15). 

The coefficient cz in (6.15) is an increasing function of P,/R} where P, is the 
number introduced into the calculation by our use of Porncarf&’s inequality 
(5.11). Now it can be shown that, for a sphere, 


(6.16) 


We say a word about how this number is obtained. It can be shown [4] (Chap. VII, 
§ 7) that P5*, which is the infimum of 


Shifi dv 
Sé 


So 
for all square-integrable functions f with average value zero and with piecewise 
continuous first derivatives such that /,, /,;is integrable, is the smallest positive 
number yu? for which there exists a square-integrable function g, with piecewise 
first and second derivatives, such that g ; g ; is integrable, satisfying 


Veg+ wg=0 
in ©;, with zero normal derivative on the boundary. By the method of eigen- 
function expansions, it can be shown [4] that g is a product of three functions, 
one each of the three spherical coordinates, and in particular, the radial function, 
say o(r7), satisfies the spherical Bessel equation 
5 n(n-+1) 
TS A4 drat te 7? 


O=0 


with an integer, do/dy=0 on the boundary and 


R6 R6 


fordr<aan,  f r(doldr)*dr< oo. 
6 6 


From this it follows that w is the smallest positive number satisfying tan (uw F,) 
= (Rs), from which uw? R}< 20.19, and hence (6.16). 
Using (6.7) and (6.16), we get from (6.15), 


fort dvScy(d) f @dv, 
S S : 

2 (6.17) 
2 7 2a) 
3 | (20.19) (2a3—3 A?+1) 
Inequality (6.17) holds for all A>41, since we have left Rs arbitrary, but smaller 
than R. Thus, in particular, it will hold for the smallest value of c,(A), A>1. 
The minimum of c,(A), A>1, occurs at about A=1.4 and gives an inequality 
for the sphere of the form (4.3) for all "w satisfying (4.2) whose second derivatives 
are continuous in G, with 


¢4(A) = 


eat] (6.18) 
to the best integer. 
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Now we remove the restriction that the first and second derivatives of 4, 
be continuous in G, maintaining the side conditions of the main case, (1.5)3. 

Let G, be a sphere concentric with © and of radius R,<R. Then the first 
and second derivatives of u, are continuous in @,. Thus defining 


o 


= = =: 1 nl 

O=0j0;, O7=T wo, 40, V, =a, (6.19) 

o o a Go : 
Se Soe 


and defining 


U; a U, — Oj x; 6 (6.20) 


we see that the corresponding strain and rotation measures are given by 


Fd; = dij, jj = Vij — Wi; (6.21) 
and hence, using (6.19) and (6.21), 


f*o,,dv=0. (6.22) 
Ss 


Thus from our foregoing analysis we get 


ftotdv<K” [@dv<K" f @adv, (6.23) 
Soa So S 


and hence, using (6.21), (6.22) and (6.23), we obtain 
aie a 2d diss 2dvu+V, 0° w* dv 
Je v ae Dare, dee es |) ody + Opes hs 
<K" [@dv+Voq?+ f wdv. 
S Ke S—GSog 


(6.24) 


But since f w,;; 4v=0, @? may be made as small as desired by choosing R, close 
fo o 


enough to Rk. Hence (6.23) implies an inequality of the form (4.3) for all a; 
satisfying (4.2) with K’’=17. In other words, 17 is an upper bound for Korn’s 
constant in the main case and, from (4.11), 20 is an upper bound for Korn’s 
constant in the second case. The upper bounds announced in (1.9), are established. 


7. An upper bound for Korn’s constant for the circle in the main case 
For the circle, using (5.3) and (5.4), we obtain 


v.17 = B[(2R/P) nn; + (1— R) 4], in C—G, Va 
and 


Vij =—4(A2—1)6,;, in Cp. (7.2) 
From (7.2) and 4° follows 
Cy pe i Oa (7.3) 


As in the case of the sphere, the minimum of y in Gy is assumed on G; and is 
given by 


Cy RS = (4) (1 — A? + A? In A2) R2, (7.4) 
Thus, using (5.12), cs becomes 


; 422 6 
oa 4—A2+22In 72 ay (7.5) 
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where P=P,/R3. It may be shown that P is independent of R,, either by simple 
dimensional considerations or by evaluation of P,; in a manner similar to that 
used for the sphere. 


Substitution of (7.1), (7.2), and (7.3) into (5.19) yields 
for dv Sey f [= (e[2) 4 — (a2) (tr d)?] dv + 
€ é 


c rd)? 2l72) A..n.n,\dv ve) 
+6. f (ltr)? + (Rip) A,,nn,)do, 


where 
A,, = — (trd)?6,,; + @6,;+2d;;d,, — dj, d;,. (7.7) 


u 


Now to our surprise and delight, A;;=0, which can be shown most easily by 
referring all quantities to a frame of reference in which d;; is diagonal. Thus, 
observing that the first integrand on the right-hand side of (7.6) is non-positive 
and using (3.3), we get 


fords | eon CSE 4) J @dv. (7.8) 
Cc € 


Since (7.8) holds for all A>1, it holds as A->oo, in which case we get 
fordvs4fdidv. (7.9) 
€ € 


The restriction that the first and second derivatives of u; be continuous in € 
may be removed in the same manner as for the sphere. 

Thus we have proven that 4 is an upper bound for Korn’s constant for the 
circle in the main case and, by (4.11), that 6 is an upper bound in the second 
case. All of the bounds announced in (1.9) are now proven. 

We close with some remarks on the extension of the above analysis for the 
circle to more general plane regions, namely, those onto which © is mapped 
conformally by a function ¢(z)=&+77, z=x-+17y, which is regular in ©. If we 
denote such a region by ¢(€) and express y, defined in the x, y plane by (5.3) 
and (5.4), as a function of € and 7, we may note the following: 

a) The boundary of £(€) is £(6). 

b) Since y, and y, are linear combinations of y, and p,,, it follows that 
y and its derivatives are both zero on the boundary of ¢(€). 

c) The Laplacian of wy transforms according to 


C'(2)|*V? y. 


d) Because of the regularity of ¢(z) in ©, |¢’(z)|? has a positive maximum and 


minimum in ©. Thus there are positive numbers, say 0 and B, bs and Bs such 


tae Pe oe bo Gn Z(G = 6) 
— B,<V?y<—b; in C(&,). 

Thus the function p defined in the &, 7 plane satisfies conditions 1°, 2°, 3°, 
and 5°, but 4° must be replaced by d°. Using p(&, 7), then, we may carry through 
an analysis for €(€) similar to that in Section 5, except that equations and in- 
equalities which follow from 4° must be replaced by equations and inequalities 


implied by @°. 


Aes 
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On the Global Asymptotic Behavior of Nonlinear 
Systems of Differential Equations 


Jee) EE VIN 


Communicated by L. CESARI 


1, Introduction 
We consider the real system 


(1.1) &=Pplt,x,y), y= Q(t,%,y) (=dldt), 


where %¥=(%1,..., %,), V=(¥1, -++» Ym)» Our primary purpose is to find conditions 
under which all solutions of (1.1) tend to zero as t-+co. When this situation 
prevails we shall, for brevity, speak of (1.1) as being ‘“‘globally asymptotically 
stable” even though this terminology usually denotes a stronger assertion about 
the solutions of (1.1). (Note, for example, that x =0, y=O is not assumed to 
be a solution of (1.1).) The conditions that imply global asymptotic stability 
will arise quite naturally from an investigation of the more general problem of 
the asymptotic behavior (as ¢->oc) of all solutions of (1.1). 

The splitting of (1.1) into two parts anticipates the key hypothesis. It is 
assumed that there exists an “‘energy’”’ function E(é, x, vy) which is positive 
definite on R,.,, (see §2 for terminology) and whose total derivative with 
respect to (1.1), denoted by E’(#, x, vy) where 


m 


(1.2) E(t, x,y) = Et, x,y) +2, Ex lt, %, 9) Bilt % y) +2 Ey lt, % 9) Gt % 9), 
iw 


jt 
satisfies 


(1.3) E'(i, x,y) S— Wy) +e (i) +e0(t) El xy), 


where W(y) is positive definite on R,, and the e,(¢)E€L,(0, co). The e,(¢) are 
introduced in (1.3) in order to permit the existence of integrable forcing or 
perturbation terms in (1.1). It is seen from (1.3), in the case e,(t) =0 (k=1, 2), 
that E(t, x, y) is (in general) negative definite only on a subspace of R,4,,. 

The existence of an energy function which is positive definite on the entire 
phase space and whose total derivative is negative definite (or almost negative 
definite, in the sense of (1.3)) only on a subspace is precisely the situation that 
occurs in several applications. An important example is the Li¢énard equation 


(1.4) x + h(t, x, x)x + f(x) =e(t), 
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where hy is a positive damping coefficient; fis a spring, that is x f(x) >0 (++-0% 
and e an integrable forcing term. If one writes (1.4) as the equivalent systema) 


(1.5) s=y, y=—hit,x,y)y—t(x) +e : 


and introduces the energy furiction 


\ 


E(x, 9) =x) +4 (ee = S76) 48,), 


then E'(t, x, ¥) = — h(t, x, y) y? + ye(t). 


If h=k>O, where k is a constant, then it is easily shown that 


E'(t, x,y) <— Wy) +|e(| + 2/e| £9), 
where W(y)=k y?. 
In terms of the autonomous case of (1.1), 


(1.6) 4= Pr PSG), 


the main lines of the argument are as follows. Under appropriate hypotheses, | 
including local existence but not uniqueness of solutions of (4.6), it is shown 
for any %, Yo that any solution x (é), y(t) of (1.6) satisfying x(0)=%9, ¥(0)=¥o: 
exists on 0<t< 00; is bounded there, the bound depending only on %», V9; and 
y(t) >0 as too. With x(t) thought of as a curve in R,, it is then shown that 
the cluster points of x(¢) as too are contained in the set Q, defined by 


(1.7) 2, = {| q(%, 0) =O}. 


(If, in a particular case, Q, contains only the origin, then it necessarily follows 
that x(t)—>0 as too. This is easily seen to be the situation for the autonomous 
case of (1.5), where q(x, 0) =—/f(x) and /(x)=0 implies %=0.) Finally, with 
the aid of the system 


(1.8) ~=f(9,0), 


it is shown that the cluster points of x(¢) are contained in a set Q, which is 
independent of %, yo and is (in general) a proper subset of 2,. In particular, 
any point of 2, which is ever mapped out of Q, by a solution of (1.8) cannot 
be a cluster point of x(t). Alternatively, the cluster points of x(¢) are contained 
in the set of points, 2, through which there passes a solution of 


(1.9) P=P(y,0), O=4(9,0). 


Thus, taken together with the other hypotheses, a sufficient condition for the 
global asymptotic stability of (1.6) is that p(t) =0 be the only solution of (1.9). 
Besides (1.4), other interesting examples for which the present methods apply 
occur in the stability theory of heterogeneous nuclear reactors, see [4]. Once 
the appropriate energy functions are shown to exist, Corollaries 1 and 2 below 
contain the results of [4] as very special cases. When restricted to the systems 
of [4], Corollaries 1 and 2 also improve the earlier results in that less smoothness 
hypothesis is required, e.g., the Lipschitz condition imposed on the [(x)-of (4.4) 
in [4] may be replaced by continuity. In fact, it may easily be shown that 
Corollary 1 also permits the f(x) term in (1.4) to depend suitably on ¢. 
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It should be observed that the crucial hypotheses of Corollaries 1 and 2, 
0=, and 0= respectively, can often be verified easily by inspection. The 
systems of [4], as well as many other easily constructed examples, are of this type. 

Lemma 3.1 is a boundedness result for nonlinear systems of differential 
equations. When restricted to (4.4), it is an improvement of a theorem due to 
ANTOSIEWICZ [1]. If n=0, 1.¢., if (1.4) is replaced by y=q(t, y), then Lemmas 3.1 
and 3.2 supplement Theorem 1 of ANTOSIEWICz [2]. 

The technique of analyzing the stability properties of differential equations 
by means of energy functions is often called LyapuNov’s second method. A 
recent survey of this method is [3]; see § 4 of [3], where there is some discussion 
of global asymptotic behavior. 

2. Summary 

The following notation is used: I,={t|tj<t<oo}, where 420. I=Ip. 
R,, k-dimensional Euclidean space. |«|=)'|x,| for x€ R,. d, the metric induced 
on R, by the norm | |. K, K,, K(), finite positive real numbers (which may 
vary from line to line). SF = fu] x Ry, [2 Sais OC), Cue ea lsIp ia), te, sets. Ol 
(vector) functions that are, respectively, continuous, continuously differentiable 
with respect to each variable, satisfy a Lipschitz condition (on the appropriate 
space). L,(Z), the set of Lebesgue integrable functions on J. A scalar function 
F(x) is said to be positive definite on R, if F(x) >0 for all «+0 and F(0)Z0 
F(t, x) is said to be positive definite on FR, if there exists an (x) which is positive 
definite on R, and F(t, x) =F (x). 

It is assumed that # and g may be written as 

P=Per Pp. I= Ga + Up; 


where #, and gz may be thought of as integrable forcing or perturbation terms, 
and where the following smoothness and growth conditions are satisfied: 


(2.1) p9EC(EXRy 4) 

For every K: 

(2-2) p is bounded on IX Sey; 

(2:3) q is bounded on I x Sk,.,,; 

(2.4c) there exists py (x) €Lip(S%) such that 
tim pelts 29) = Balt), 


where the limit is taken in t and y together, and is assumed to be uniform 
with respect to x on Sk; 


(2.48) there exists a scalar function p,(t; K) €L,(L) such that 


(2.5a) there exists g,(x)€C(R,,) such that 


lim da(t, *, ¥) = 44(%), 
t— oo, y> 0 
where the limit is taken in t and y together, and is assumed to be uniform 
with respect to x on si 


5* 
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(2.58) there exists a scalar function q (t; K)€L,(L) such that 
last, x 1S a(t K) (4%, 9) TCS) Fan 


It should be noted that in the autonomous case p, =0, gg =0 so that (2.48, 2.58) 

are trivially satisfied; (2.1) implies (2.2, 2.3, 2.5); and PE C’(R,4 4) Is*suiieient 

for (2.4). 

Concerning the energy function E(t, x, y) and the associated E’(t, x, v) defined 

by (4.2) it is assumed that: 

(2.6) There exists E(t, x, y)€C'(IXRy4m) and a positive definite E,(x, Vv) EC (Ry+m) 
such that E,(x, y)SE(t, x, y) and E,(x, y)—0e as || + |y|—oo. Further- 
more, E(t, x, y) 1s bounded on Ix St, for every K and E(t, 0, 0) =0. 


(2.7) There exist e,(t) €L,(L) (k=1, 2) such that 
E'(t, x,y) Sa (t) + a(t) El, *, 9). 


(2.8) There exist e,(t)€L,(1) (k=1, 2) and a positive definite W(y)CC(R,,) such 
that W(0)=0 and 


E'(t, x,y) <— Wy) +er(t) + ee(t) Ele x 9)- 


Although (2.7) is implied by (2.8), it is listed separately, since it is sufficient 
for Lemma 3.1. We shall assume, without loss of generality, that the e,(¢) of 
(2.7, 2.8) are non-negative. 

Let y, %o Yoel XRyag BY SHLG te: Yo We Meanaany -soluuen 
x(t), y(t) of (1.1) that exists on some subinterval of J having ¢) as left-hand 
end point and that satisfies x(t) —=%», V(ty)= Vo. It is seen from (2.4) that for 
any given (tf), %, Yo) an Y exists but is not necessarily unique. w€R,, is called 
a cluster point of /, for an / existing on J,,, if there exists a sequence {t,} 
with ¢,>0o as k->oo such that x(4,)—q@ as k—oo. More properly, we might 
call such an w€R,, a cluster point of x(¢) rather than of Y. The definition anti- 
cipates that, under the preceding conditions, for any ¥, y(t)—>0 as too, as 
will be shown. Thus, (x, y)=(q, 0) is really a cluster point of Y. The set of 
cluster points of Y is denoted by /*. 


Central to our considerations is the subset Q, of R,, defined by 
(2.9) G24 Gs (ae) Oh 
From the continuity of g, it follows that Q, is closed. 


Theorem 1. Let (2.1, 2.2, 2.3, 2.5, 2.6, 2.8) be satisfied, and let (fo5.%os00/o) 
CLXR im. Then any S= St, ty, Xo, Yo) exists on I,, and (x(t), y(t)) € S¥.4y for 
some K=K (ty, %9, Yo). Furthermore, y(t)>0 as t>co and S*CQ,. 

In Lemma 3.4 it is shown that 0€2,. The following result on global asymp- 
totic stability is an immediate consequence of Theorem 14. 

Corollary 1. I/ im addition to the hypothesis of Theorem it is further assumed 
that 0O=2,, then it also follows that x(t)>0 as t>o. 

The conclusion £* CQ, of Theorem 1 can be considerably improved if con- 
dition (2.4) is also assumed. For this purpose, we consider the system 


(2.10) p = pa(@): 
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Let p(t, €) denote the solution of (2.10) satisfying p(0, )=é. It is well known, 
note the Lipschitz condition of (2.4«), that the mapping /, of R,, into R,, defined 
by 


(2.11) Mé=9(t8)  (E€R,) 

is locally a homeomorphism. The set 2, <Q, is defined by 

(212) Q, = {x|*%EM2, and, forsome t=i(x)>0, M,xGQ,}. 
Define 

(2.13) Q=Q,—Q,. 


It is clear from the preceding definitions that 2 may be equivalently described 
as the set of points in R,, through which there passes a solution of 
(2.14) P=bal~), 0= Gal). 


Theorem 2. Let the hypothesis of Theorem 1 as well as (2.4) be SBR eee Then 
it also follows that S* CQ. 


In Lemma 4.3 it is shown that 0€ 2. The following result on global asymp- 
totic stability is an immediate consequence of Theorem 2. 


Corollary 2. Ij in addition to the hypothesis of Theorem 2 it ts further assumed 
that O=Q, then it also follows that x(t)—>0 as t>ov. 


3. Proof of Theorem 1 


Lemmas 3.1, 3.2, and 3.3 establish Theorem 1. Lemma 3.4 shows that 0€2, 
as asserted in § 2. 


Lemma 3.1. Let (2.1, 2.6, 2.7) be satisfied, and let (ty, %, Yo) CL XRy4m. Then 
any S=Lt, ty, Xp, Vo) exists on I,, and 
(3.4) Ie@| SK, |ly@[/SkK (¢€l,) 
jor some K=K (ty, Xo, Yo)- 
Proof. By (2.1) there is a 4)=4,(Y%)>¢% such that # exists on 4StSh. 
Define 
ny aa 
th, 


E(t)=E(t,x(t),v()), E,(j= 
63.2) (?) (t, x (4), ¥(¢)) 


(x(¢ 
E' (t) = E' (t,x), vy) ee Sis 
Then (1.2, 2.7) imply 
E()=E'WSat+eQZO (StS), 


and thus 


E(t) SE t) + Fel Jar fete Laide Whii=s 74) 
A well known inequality now yields 


63) HOSE) S{Et i) + Fa 1) dt} exp| fea(0) dt] (tyS¢S4). 


to 
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Upon noting that the right-hand side of (3.3) is finite and independent of 4, 
the conclusion follows readily from (2.6, 2.7, 3.3) and a routine continuation 
argument. 

Lemma 3.2. Let (2.1, 2.3, 2.6, 2.3) be satisfied, and let F be as in Lemma 3.1. 


. Then 
(3.4) lim y (¢) = 0. 


t— co 


Proof. Define E(t), E’(t) on I, by (3.2). Then by (1.2, 2.8) 
EB (i) =E'() S— Wy) ta +a FY, 


from which 
3.5) OSE(t))S E(t) —fW(y(t)) dt PAR NOES AEE ICS: 


It follows readily from (2.6, 2.8, 3.1, 3.5) that W(y(t)) €L,(J,,). 

Suppose (3.4) is not true. From (1.1, 2.3, 3.1) one has |y(#)|SA(¢CJ,) for 
some K. Hence, there exists an ¢,>0 and a sequence of intervals {(t,,, t,+A)}, 
where t,, >0o as m->co and A> 0, such that 
(3.6) Wet By AG, ete Ae Ane 
But from (2.8, 3.1, 3.6) one has for some ¢,>0 

W (y(t) = & @, Sts +7; w= 1,2)...)5 


which contradicts W(¥(#)) €L,(J,,) and completes the proof. 


Lemma 3.3. Let (2.1, 2.2, 2.3, 2.5, 2.6, 2.8) be satisfied, and let S be as im 
Lemma 3.1. Then 


(3.7) ae nO: 
Proof. Let A>0 andi#¢€J,,. Then from (1.1) one has 
t+Aa mn tra 
y+ A) — 9) = JS {galt ¥(t), ¥(t)) — Gala (t)) fae + J Ga (x(x) d+ 
t 
tA 


+ f 9p (, ¥(2), (2) dz, 


which together with (2.5, 3.1, 3.4) easily implies 
t+a 


(3.8) lim IR Clea cha 10 
OOF 
We wish to show that 
(3.9) jim q(x (4) =0. 


Suppose (3.9) is false. Then it is false for some component, say 9, i(x(). Then 
there exists an ¢,>0 and a sequence {i,}, where t,—>co as h->0o, such that 
| Ga, « (¥ (44)) |= & (R=1, 2,...). We assume that Iu, i (X (tx) = &, a similar argument 
holds for ¢,, ;(«(&)) << —e,. Uniform continuity, which one has from (2:5 005.34) 
implies the existence of a 6>0 and independent of k such that 


(3.10) Gi(6)2e% (1E— (4) SO; k=1,2,...). 


? 
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However, by (1.1, 2.2, 3.1, 3.10) it follows that there exists a 4,>0 and inde- 
pendent of & such that 


(3.14) (A) 2ee —~GStSt,+A5k=—1,2,...). 


Clearly (3.11) contradicts (3.8) for A=A,. Hence, (3.9) is true, and this easily 
implies (3.7). 

Lemma 3.4. Let (2.1, 2.2, 2.3, 2.5, 2.6, 2.8) be satisfied. Then q,(0)=O, 07, 
equivalently, OC Q,. 


Proof. Let ¢,>0 and ¢,>0. Then by the continuity of 9, it suffices to 
show that there exists a CR, such that |&|<e, and |g, (&)|<e,. By (2.6, 3.3) 
tj may be chosen sufficiently large so that for any A= A(t, ty, 0, 0) it follows 
that |x (2)| <¢, for ¢€I,. From (3.9) one has, however, the existence of a 42% 
such that |¢,(%(4))|<e. The result follows immediately on choosing €= x(t). 


4. Proof of Theorem 2 
Lemmas 4.1 and 4.2 establish Theorem 2. Lemma 4.3 shows that 0€2 as 
asserted in § 2. 
Lemma 4.1. Let (2.1, 2.3, 2.4, 2.6, 2.8) be satisfied, and let S be as in Lemma 3.1. 
Let e>0 and A>0. Then there exists ann(é, A, ty, Xp, Vo) > O and aT, (e, A, S) CI, 
such that t, €I7,, CR, and satisfying 


(4.1) |~(4) — El <9 
implies 
(4.2) Ix) —M,é|<e (4StS4+2). 


Proof. By (3.1) there is a K,>0 such that |~(¢)|, |y(¢)| SA, for ¢¢J,,. With- 
out loss of generality we assume that e< K,. Let K, denote the Lipschitz con- 
stant in S?™: of (2.4a) for £,. Choose y>0 so that 


(4.3) nS teexp[— KA]. 
By (2.4, 3.1, 3.4) we can choose J, €/,, such that 


[Pa (t, *@), 9) — ba(*)| SS exp [— KzA] (#E La), 


= 6a. 
(4.4) BS 
J p(t; Ky) dt S % eexp[— K, A]. 
is 


Let ¢,€ 7, and let & satisfy (4.1). Using (2.11), define p(?) by p(t) =M_,, 
&—qg(t—t, &). As #, is continuous there exists a ¢,, where t<tpSt,+A, such 
that | y(é)| S24, for 4S¢St,. Since (2.10) is autonomous, one has from (1.1, 
2.10) 


x(t) — p(t) = %(4) — E+ Jipa(t, x(t), y(t) — pa (*(t))$ae + 


+1 Pax) Alor + fone) yO) de StSH), 


12 Jett: Levin: 
which together with (2.4, 4.3, 4.4) yields 
t 
v(t) — (| <deexp[—KiA] + Kas |x() —e@)lde  (HStS4). 
ty 
A well known inequality now implies 
(4.5) |«(t) — p()| See exp[Ka(¢—4—4)] (StS), 
* which in turn implies 
(4.6) |p| SK, +heexp[Kelt—h—a] (4StS4). 


The conclusion follows readily from (4.5, 4.6). 

Lemma 4.2. Let (2.1—2.6, 2.8) be satisfied, and let F be as in Lemma 3.1. 
Then S*oQ. 

Proof. Suppose there exists an w€ Y* but a¢Q. Then by (2.13, 3.7) one 
hasw€S*NQ,. Since w€ F*, there exists a sequence {f,}, with t, +00 as R00 
such that 


(4.7) jim £(t,) =o. 
Since w€Q,, there exists (see (2.12)) a A>O such that M,@ ¢Q,. Hence, it 


follows, as Q, is closed, that d(M,w, 2,)=e@>0. The homeomorphic nature of 
M, then implies that there exists a v, where O0<»<¥g, such that 


(4.8) d(M,0,2)=49  (l-o| <»). 


Choose N sufficiently large so that by (4.7), (3.7), and Lemma 4.1 (with 
£=x(ty)) one has respectively 


(4.9) | (ty) —@| <z, 
(4.10) d (x(t), Q) <4, tel), 
(4.11) |x(t) —M_y»x(y)|<t (yStSty+A). 


From the triangle inequality 


d (My * (ty), 21) S| My x (ty) — % (ty + A)| + € (x (by + 2), 21) 


(which is easily shown to be valid even though one ‘‘point’’ is Q,), (4.8, 4.9), 
and (4.10, 4.11) at ¢=ty+A one has }9<4y7-++4y», which contradicts y<$o and 
completes the proof. 

Lemma 4.3. Let (2.1—2.6, 2.8) be satisfied. Then 0€Q. 

Proof. By (2.12, 2.13) and Lemma 3.4 it is sufficient to show that (t, 0)=0 
for ¢€I. This is, as (2.10) is autonomous and 4, satisfies a Lipschitz condition, 
equivalent to £,(0)=0. Hence, it suffices to show that p(t, 0)=0 on 0OStSA 
for some A>0. Let A>0 and y>0 be sufficiently small so that by (2.4«) v(t, &) 
exists on OStSA, |E| <>. 

Let ¢;>0 and e,>0, with e,<». By continuity it suffices to show that there 
exist a €€ R,, such that |é|<« and | y(t, €)|<e, on 0StSA. By (2.6, 3.3) ¢) may 
be chosen sufficiently large so that for any a SAb st. 0, 0) it follows that 
| x (t)| <min(¢,4é) for tEl,. From (4.2) one has the existence of a t,=t) such 
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that 


|x (t) — p(t —4, «(4))| << Fe, (StS4,4+A), 
which implies 

| p(t —4, *(4))| < Be (4<tS44+4), 
and, hence ‘ 

| @ (4, #(4))| < Fee (OStSA 
The result follows immediately on choosing £= x (t,). 


5. Supplementary Remarks 
While (5.1, 5.2) below are essentially implicit in the preceding sections, the 
following modification of the proofs yields a little more information. This modi- 
fication has been pointed out to us by Y. SipuyAa. Lemmas 3.1 and 3.2 are estab- 
lished as above. Then, let@w€-Y* for some S= S(t, to, Xp, Yo), and let x (t,) > 
as k->oo, where ¢, >0co as k->oo. By an argument similar to that of Lemma 4.1 
it may be shown that p(t, w) exists on J and that 


(5.1) Jim x(t-+4) = ol.) 


uniformly with respect to ¢ on any finite ¢ interval OS¢< 7. Then by (5.1) 
and the same calculation used to establish (3.8) one can easily show that 


Jala ro) aa = 0 (t€ I) 
and, hence, that 
(5.2) qa (p(t, )) = 0 (¢€1). 
Clearly (5.2) implies that Y* CQ. 

The method of the above paragraph has the advantage of showing that (2.2) 
is not a necessary hypothesis for Theorem 2. However, it appears that (2.2) is 
necessary for Theorem 1 unless, of course, (2.4) is assumed, in which case Theorem 1 
is a trivial consequence of Theorem 2. It may also be noted that (2.3) cannot be 
entirely dispensed with, as has been shania in Remark 2.2 of [4]. 

As alluded to in the first paragraph of §1, the conclusion of Corollaries 1 
and 2 is not that the origin is asymptotically stable in the usual sense, see [3, 
p- 142], but merely that all solutions of (1.1) tend to the origin as too. For 
the origin also to be stable, it is sufficient to add 


(5.3) $(t,0,0)=0, g(t,0,0)=0, e&(¢)=e,(%) =0 (tET) 


to the hypothesis of Lemma 3.1, as is evident from (3.3). Invoking (5.3) rules 
out, of course, many of the problems that are included in the present considerations, 
e.g. (1.4) with e(t)==0. For the origin also to be uniformly stable, it is sufficient 
to assume further that 
lim E£(¢,x%,y)=0O uniformly on ?¢€J. 
|*|+|y|>0 

Added in Proof. Recently J. P. La SaLeE in the Proc. N. A. S. 46, 363—365 (1960), 
has announced results concerning autonomous systems which overlap considerably 
those of the present paper when the latter are restricted to (1.6). 


The work reported in this paper was performed by Lincoln Laboratory, a center 
for research operated by Massachusetts Institute of Technology with the joint support 
of the U.S.Army, Navy, and Air Force. 
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On A pproximate Solutions 
of Non-linear Hyperbolic Partial Differential Equations 


ROBERT H. MOORE 


Communicated by L. CESARI 


§ 1. Introduction 
Consider the hyperbolic partial differential equation 


(4.4) Upy =] (%, Y, U, Uy, Uy) 
given together with the characteristic initial condition 


(1.2) u(x,0)=o(x), 4(0,y)=Tt(y) (where o(0) = 7(0)) 


to be satisfied by u(x, y) on the rectangle #: OS xSa, OS ySb. 
In comparison, consider also the ordinary differential equation 


(1.3) u! = f(x, u) 
with the initial condition 
(1.4) u(0) =o (= constant) 


to be satisfied by u(x) on the interval 4%: OS wa. 

These two problems will be denoted respectively by 2-HP and 1-Ord. 

Many authors have realized the analogy between these initial value problems 
in two and one independent variables. The object of this paper is to develop 
a proof of the existence (under suitable conditions) of a solution of 2-HP which 
is analogous in some respects to the existence proof given by PEANO [9] for 1-Ord. 
The important part of this development is that a certain broad class of “‘approxi- 
mate solution functions’ are indeed approximations to a solution in the sense 
that, under a natural ordering, they converge uniformly to a solution. The 
actual conclusion of the existence of a solution under the given conditions is 
less important, for this is already known (cf. HARTMAN & WINTNER [5], P. LEEHEY 
[7] and J. B. Draz [3)). 

A subsidiary aim, for the purposes of the present paper, will be to extract 
the essential notions in the development of the proof in such a manner as to 
suggest how a proof might be developed for analogous equations of higher order 
in correspondingly more independent variables. 

In a subsequent paper, an analog of the Runge-Kutta procedure for numerical- 
ly obtaining approximate solutions of 1-Ord. will be developed for 2-HP. The 
procedure derived there is similar both in application and order of accuracy to 


76 Rosert H. Moore: 


that given by Kurta [6]. Its validity depends on the proof in the present paper; 
indeed, the numerical process defines an approximate solution which is included 
in the class of approximate solutions investigated here. 

In viewing 2-HP as an analog of 1-Ord., one may consider generalizing the 
Euler-Cauchy polygon method of finding approximate solutions of 4-Ord. In 
this method, one subdivides the interval ¥% at the points 0=&5, Ea eee 
into subintervals 4&,=€,,,—&; and defines the function v(x) on 4% by 


(1.5) v(t) = + DC, AE, +C(x—&), where C,=7(E.,0(,)), 


for £,<*<&,,,. This defines v(x) as a continuous polygon (7.e. piecewise linear 
function) through (0, 0) with vertices where x=€;. ZWIRNER [10] gave a gene- 
ralization of this for the special case of 2-HP when u,,—/(x,y,u). This he 
did by subdividing # into subrectangles in the usual manner by a mesh of lines, 
and defining over each subrectangle a bilinear function (hyperbolic paraboloid) 
in such a manner that these bilinear functions could be pieced together to obtain 
a continuous function approximating a solution on all of #. Analytically, such 
an approximation can be formulated by choosing division lines x=€; and y=n; 


in, letting AE,=£,,, —§, Anj=njia — and G(x) =o (&) + CEW—F) (yg) 


x (ny +a)— (i) Ag, 
f earns 
for &Sx<€,41, t(y)=t (yn) + 2 — ™ (y —n) for 47S y<nj41 and finally 


setting Af 
4—1 j—1 
v(x, ¥) =a(x) +7(y) —o(0) + 2d 2 Cup 4b, Ang t 
(1.6) ie 


j—-1 
=t 20 Ane («—&) + 2 Cy; ME, (v— nj) + Gj (*— &) (y— nays 
= a=0 

where C,, is a constant approximately equal to / (&, mg, v(&,, np)). If exact 
equality were required here for C,,, then the analog with the polygon method 
for 1-Ord. would be precise; but in fact this greater degree of generality is 
permitted, the only requirement being that as the mesh is made finer the value 
C.,, Should tend to f (E.. np, v(&,, Ng)). Actually ZWIRNER gives specific restrictions 
on how C,, may deviate from /, but the requirement mentioned is satisfied 
because of these restrictions. 


Both (1.5) and (1.6) can be motivated by putting the corresponding differential 
equations into integral form and then carrying out an approximate integration 
with a piecewise constant integrand. Putting (1.1) and (1.2) into integral form, 
one finds 


x 


u(x, y) =o(x) +7(y) —o(0) +f 


Fils, t,u(s,t), P(s, t), g(s, 2) ds dt, 


(ee Biles) — 61 (x) + Ft (x t, u(x,t), P(x, t), q(x, 4) dt, 


g(x,y) =r'(y) +51 (s,¥,4(s, 9), P(s, 9), 4(s, 9) ds. 


These form the starting point of the analog of the Euler-Cauchy method developed 
recently by Diaz [3]. The integrals are replaced by summations quite similar 
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to those in (1.5) and (1.6). In particular, (1.6) yields an approximate solution 
of (1.1) by taking (cf. p. 368 of Draz’s paper) 


Cup =F Exe nes v(x 7), Voate 1) = Ubu mp) 2 (Eu, ptr) — 9 (Ea. 8) 
p= TEx ne, ¥ (Exp) pani aeeeee ) 


This choice of C,, yields approximating functions for (1.1) quite analogous to 
the Euler-Cauchy polygons for (1.3). Again the approximating surface is made 
up of segments of hyperbolic paraboloids. 

By including wu, and uw, among the arguments of /, one encounters an im- 
portant difficulty. When in the second and third of equations (1.7) the inte- 
grand f is replaced by a piecewise constant function C, the resulting approxi- 
mations to f and g are not continuous. For the approximations to #, for example, 
there is no integration with respect to x to smooth out the discontinuities of C 
on lines x=€,. This difficulty was encountered by Diaz. To overcome it, he 
derived a lemma (stated in §3 below) which extracts the needed properties. 
On the strength of its conclusion one may obtain the needed convergence of 
the approximations of the derivatives by applying a theorem of ARZELA which 
deals with the convergence of certain not necessarily continuous functions to 
a continuous limit function. J. CONLAN [2] gives a proof of ARZELA’s theorem 
in an appendix to a paper in which he extends D1Az’s method to the Cauchy 
and Mixed Boundary-Value Problems. 

In the development below for 2-HP, the domain space of independent variables 
will be divided as usual into subrectangles. Now it may be observed in both 
the Euler-Cauchy polygons for 1-Ord. and their analogs for 2-HP that the most 
important values of the approximate solutions are those on the nodes of the 
meshes. The values at intermediate points are unimportant. For example, for 
{-Ord. the polygonal approximation (1.5) could be replaced by a step function 
agreeing with the polygon at the left endpoint of each subinterval of J. The 
sequence of step functions so obtained still converges uniformly to the solution. 
Therefore, so far as existence proofs and numerical processes are concerned, the 
values at non-node points as given by a continuous polygon are of no concern. 
In this presentation we have thus omitted any consideration at all of the inter- 
mediate values when defining the functions approximating the solutions. Upon 
defining these functions in the usual stepwise manner beginning with the initial 
conditions, one may observe further that they satisfy certain ““growth conditions” 
which appear in the form of summations resulting from approximate integration 
of the integral equations equivalent to the given differential equation. These 
growth properties lead, via the inequality in D1az’s lemma, to the satisfying of 
the ‘‘smallness of discontinuity’? requirement of ARZELA’s theorem, which in 
turn assures the existence of the desired subsequence of approximate solutions 
converging uniformly to a solution. 

In accordance with the above, in §2 ARrzELA’s theorem is given in a gener- 
alized form for functions whose domains are the node sets of convergent meshes 
in Euclidean n-space and whose ranges are in a complete metric space. In §3 
the concept of a “‘vine”’ is introduced, a vine being a function defined on mesh 
nodes whose image at a given node depends only on the images at preceding 
nodes. The basic result of the section (Theorem 2) shows that for sequences of vines 
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defined on convergent meshes and satisfying certain growth conditions there are 
subsequences converging uniformly to continuous functions. In § 4 it is shown 
through two lemmas that sequences of vines of the type considered in Theorem 2 
and appropriately related to the problem 2-HP do exist, and that the limit 
functions so determined yield solutions to the problem. These results are sum- 
marized in Theorem 3, which is our principal result. The final theorem, given 
in §5, shows that every continuously differentiable solution of 2-HP is the 
limit of a sequence of approximating functions of the type under consideration. 
This answers affirmatively, for u,,=f(x, y, 4, ~, 7), a question suggested by the 
results of ZWIRNER [10] concerning the equation u,,=/(%, y, 4). 


§ 2. Functions on Mesh Nodes 
Throughout the following 7 will have the range 1, ..., ”. 
In n-dimensional real space let 


45% denote thé point (7), 47. 55,.4 5 
n 


2. |x| = 21 4"1; 


3. # denote the cell 0S x'<a' for some fixed values a’; 
4. {@"} be a sequence of meshes on 2 formed by the hyperplanes 


x= EE” (k=0,1,...,m'”) (where &”=0 and £%,—a’) 
chosen such that o, = fue — ake is positive; 
5. WM” denote the set of nodes of .@’, that is W’= {( ae ER, the es 


6. & ky...k, denote the node ( of sie Red ae ENS 


weey 


Thus, if lim ||.” || =0, then the diameter of the largest subcell of # ap- 
proaches zero as y->oo, and the points in UW” are dense in &. 


Definition. A sequence of meshes {.@"} such that lim||./”||=0 will be 
called convergent. : 

Let ¥ be a complete metric space with distance function g. For a given 
sequence {.@} of meshes in 2 let {V"} be a corresponding sequence of functions 
VN sy. 

Definition. A sequence {V’} of functions V’: W’—>W will be called Arzela 
quasi-continuous, or an Arzela sequence, if and only if given ¢>0 there exist 
6, and N, such that y > N, and 


| Sigs a &, ke 


bt 
Rye Ry 


<0) 
imply 
O(V" (Ei, tn)» V" (Ek, tg) <e 


Note that this property of sequences is less restrictive than equi-continuity 
in that, even when an analogous definition is made for functions ¥+YV, the 
functions need not be continnous. It does require, however, that any cont 
tinuities become small as y->oe. 


t 
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Definition. The functions V” will be said to converge uniformly to a function 
V of 2 into Y if given ¢ > 0 there is an N, such that » > N, implies @ (Ve( Eon) 
V (&%,...m)) <e for each & 4, €4”. 

Although the functions V” are defined only on -””’, there is a simple extension 
to all of @: If &<x'< &?., for each 7 (this includes the possibility x*=a’), 
then let Veni) ( Eo), Note that the uniform convergence (in the usual 
sense) of V*” to V implies the uniform convergence (in the above sense) of V” 
to V. This extension is useful in the proof of the generalization of ARZELA’s 
theorem (Theorem 1 below), but it is not otherwise needed in the sequel. 

In case # is simply a rectangle in E* and Y is a linear space, another way 
to extend the definition of ¥” to all of % is bilinearly over each subrectangle- 
(cf. Diaz [3]). 

Theorem 1(ARZELA). Let & be a cell in E”, and let {M’} be a convergent 
sequence of meshes in #. Let N” denote the nodes of M’. Let V be a complete 
metric space and {V"\ an Arzela quasi-continuous sequence of functions V’: NV” > 
such that the ranges of all the V” lie in a compact subset of V. 

Then there is a subsequence {V\ of {V"} and a continuous function V: 2 >W 
such that Vi converges uniformly to V as 7+. 

This theorem is easily proved by first employing the extended functions V*” 
above mapping 2 into WY. For the sequence {V *’} the proof follows in a fashion 
almost identical with that for a sequence of equicontinuous functions, namely, 
via the diagonal process employed in connection with the countable dense subset 
of & given by the separability. Having found in this manner a subsequence 
{V+ converging to a function V, one employs the corresponding subsequence 
{V"} and the same function V to prove the theorem. (See ConLaN [2], Appendix 2, 
for a complete proof for the case of functions defined on the whole of any compact 
separable metric space 2.) 

Note that this theorem provides a continuous function on all of 2 as the 
limit of a sequence of functions each defined on discrete point sets. It serves 
thus as a bridge from functions with discrete domain to functions with continuous 
domain. In numerical work, it is the former which are actually employed as 
approximate solutions of differential equations. In succeeding sections it is 
functions of this type which are studied. At the end, Theorem 1 is used to 
provide a continuous function which is a solution of 2-HP. 


» 


§ 3. Vines 

In this section we consider specific functions V” defined on 4” which are 
suggested by the observations in the Introduction. 

Let ” be partially ordered in the natural manner: &; |», S &...4, provided 
hk; Sk; (¢=1, ..., n). The symbol < will mean both S and = hold. 

Definition. Let V” be as above, a mapping of WY” into Y. V” will be called 
a vine if the value V’(&, _;,,) depends on the values V"(&, 4) for Ge yw < Eh, ohn 
but is known as soon as these values V"(&%-_~) are known. 


For notational convenience let Vj), =V" (Eh, ... iy): 
The conditions on V” now to be considered are given for % in 2-dimensional 
real cartesian space R?; they are motivated by the second and third of the integral 
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equations (1.7). The statements following are thus designed to apply to the 
special case 2-HP of present interest. However, it 1s hoped that this inter- 
pretation is suggestive for higher dimensional cases. 

Suppose, for some constant L=0 and some quantities &,’ and ¢, 20, that 
the vines V” satisfy, for »>Ng and ||é,— &y7|| <6¢, the inequalities: 


1-1 

Gy: 0 (Vai, Ver) S exw +L DO (Ve, Vir) A PERO Mesos 
j=0 
k—1 


CP ote, Vin = a HL Do ae, BO. 
41=0 


R —1 
where, by convention, sums >) are zero. The vines V” will then be said to 
0 


satisfy the growth conditions G, and G,. The quantities e,,,and €);, can be thought 
of as initial oscillation bounds of V’. Such quantities, in the application to 2-HP, 
will be related to the initial conditions o and tT. 

With the preceding growth conditions as background, a lemma due to D1Az 
([3], p. 373) will now be stated. This lemma is the key to proving that the 
sequence of vines is Arzela quasi-continuous and subsequently that it converges 
to a continuous function V of # into v. 

Lemma 1. // 

(a) the interval [0, A] is subdivided into finitely many subintervals Ao, Ay, ..., 
Avent 

(b) £9, 1, +--+, & Us @ set of non-negative veal numbers; 

(c) the numbers L=0 and ¢=0 are such that the inequality 


s—1 


QSL Lye A, 


v¥=0 
holds. for s=A, sit; 


then the inequality 
SL 
S41] Ble LAI {e+L A} S ge {6 they alg} 


r=0 
HOUdS far SS=1K...8, fb. 
Remark. If, moreover, gy<e, then g,<ee"4(1+LA,) for s=1,2,...,t. 
This is an immediate consequence of the lemma. 
The last result applies immediately to vines V” satisfying the growth con- 
ditions G, and Gy. By setting [0, A] =[0, a®], A,;=A?", e= ey and g,=0 (Vx, Vz,) 
it follows from G, and the lemma that, for »>N, and l| $22 —Se-r || <dg, 


, 


o (Vii, Vir1) ES En Re’ en" (4 -++ oe ao’) ‘ 
Similarly, from G, it follows that 
o(V;", Viv) < Ely em (4 -- IE, ZA) ; 


If the sequence of meshes {.@’} on which the V” are defined is convergent 
then there is an N, such that » > N, implies LA}’<1 and L A2’<4. Thus, for 
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enn (M., NG), 0 (Vii, Vi v) 32 (Ey oe + Ely eet) . 


Suppose further that ¢,,, and ¢,, are such that given e> 0 there exist 6, and 
N, such that the conditionsy > N, and || &%,— || < dgimply 2 (egy 6" + ep eh) <e 
on all of W”. Then fromy > N,=max(N;, N,, Ng) and || &;—& 7 || <6, =min (dg, 5.) 
it follows that 
o(Ves, Vir) <3 


that is, the sequence {¥”} is an Arzela sequence, and the Arzela theorem immedi- 
ately applies. The foregoing results may be collected into the following theorem. 
Theorem 2. Let & be a cell {(x1, x*): OS x'<a‘} in 2-dimensional real space; 
let V" be a complete metric space. If 
(a) {&} is a convergent sequence of meshes on & ; 
(b) {V"} is a sequence of vines, V?: VW’ >V; 
(c) the V” satisfy the growth conditions G, and Gy; 
(d) the tmitial oscillation bounds &, and ey are such that given e>0 there 
is an N, and a 6, such that »y>N, and || b,,—€&yr||<6, imply exytey<e; 


then there is a subsequence {V"i\ and a continuous function V of & into V such 
that {V"} converges uniformly to V as 700. 


§ 4. Approximate Solutions, and their Convergence, for 2-HP 

In the application of the preceding to differential equations, satisfaction of 
the growth conditions depends on a function f=/(x,v) of 2 x¥ —R, the real 
field. In the present case this is the function / appearing in equation (1.1). In 
order to develop the relationship between vines and the problem 2-HP let: 

1. € denote the cell 0S x'<a‘ (i=1, 2) in 2-dimensional real space, and let 
the notation of §3 apply to #; 

2. V be 3-dimensional real space with points v= (v®, v!, v?) ; 


2 
5. lol= > 
j= 
4. {My} be a sequence of (infinite) meshes on Y formed by the planes 
vi =a)” (p= ORs 1 2 ee) 


chosen such that w/”,, — i” is positive; 

5. @},,,7, denote the point (a7”, w,,”, w7”); 

6. || -4-||= lim sup ||@%,41,7,41,7.41 —%s,n,r,|| When this exists; 

fa be en meh on & XV formed by -@ and -47;- together; 

8. ||. a*"|| = max (|||, || 4|) when this exists. 
As with.@’ ,the sequence of meshes {.@*”} will be called convergent if || #*’||>0 
as y-> oo. The sets {(x, v)} of ¥ XV such that sx ee Aner or, "Sv <oltyy will 
be called cells of M*”’. 

If V°”, Vv)” and V®” are real valued functions on WV”, let 


Vii = VP (Eq) = VIE EP) and Vin = (Ver, Val, Val) 
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so that V’ is a function on W” into Y. (In the following, V” will be a vine 
mapping V” into ¥.) 

For each y, let @” be a real-valued function, ©” (x, DC SIO 0) 
defined on & XV such that @’ is constant on each cell of *”. Here, as later, 
a point in #x¥W will be denotéd by (x, v) ‘or. (a1, 377.0% 07) 26 notational 
convenience dictates. 

Lastly, let Chy= C" (Gi) =© (Shr, Van). 

Lemma 2. (Existence of vines satisfying the growth conditions.) Suppose o(x*) 
and +(x?) ave continuously differentiable on the intervals (0, a*| and [O, a*], re- 
spectively, and o(0)=1(0). Suppose f is a real-valued function f(x, v)=f(x*, x*, 
v, v', v2) defined on XV which is continuous and bounded and satisfies on 
Lx the Lipschitz condition with Lipschitz constant L: 


Cae cast avai Da Cea rie es |e L{| yi — p| 4+ |v? — v|}. 


Let M&’ and M*” be convergent meshes on XK and XXxYV, respectively. Let 6” and 
C” be as above, and suppose further that 6 converges uniformly to f as y—>ov. 

Then there is a sequence of vines V” of N” into V satisfying the hypotheses of 
Theorem 2. 

Note. In a later use of the lemma, three sequences {61 (7=0, 4, 2) will be 
used to define respectively the V°’, V'’, V?” appearing in the proof below. 
See the remark preceding Theorem 3, § 4. 

Proof. Let o,=0(é;”) and ty=1(€é?’). Define o;,” and 1,” similarly. For each 
y define 


k—-1 1-1 
Vet 0p] og) Capa age 
e=0 p=0 
w1 
) ly ” 
(4.1) Va —— On. a ps Che Ae 
p=0 
k—1 
) ty D , 
Vel te ee Ce 
a=0 


= 

where by convention, )}=0. Set V’=(V°", V", V*"). The sequence {V"} will 
0 

be shown to be as desired. 

For notational simplicity, the index » will be suppressed in the remainder 
of the proof, but, for later use, given ec, will be thought of as so large, say y> No, 
that |@"(%, v) — f(x, v)| < ee. 

It is easily checked that V is well defined. Furthermore, V is a vine. 

Hypotheses (a) and (b) of Theorem 2 are thus satisfied; it will be shown that 
(c) and (d) are also satisfied. 

First, note that V is bounded. For, since / is bounded, @ is also bounded, 
say by M, while, because o, o’, t and t’ are continuous on compact sets, they 


too are bounded, say by the same M. Thus, it is clear from the defining equations 
(4.1) for V that 


[V1 <3M+M aa? [V4])SM+Ma? and lV?]}<M4+Ma. 
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Hence, / is considered only on a bounded portion of 2 x¥W, and, being con- 

tinuous on all of 2 x¥, is thus uniformly continuous in the domain of interest. 

Also, o, o’, t and 1’ are uniformly continuous on [0, a'] and [0, a?], respectively. 
Next, 


2 
(4.2) i ar Vall = (Vi — Vii, 
j=0 


in which, with k*=max(k, k’) and k,=min(k, k’), one has from (4.1), 


1 nk al 


Ve — Vii] S|eu— acl + 2i,|Cael 4a Ap, 


1 ait , , 2 
er ee er ee B 
=0 
k*¥—1 : 
74 2 \) 
Viiv— Ves) SD. | Curl De: 
a=kh, 
Because |@|< M ando is uniformly continuous, it is clear that for /=0, 1, ..., m?” 


both | Vin —Vj",| and | Vz; —V;7,| can be made less than a given eg by taking | & — &}| 
less than some dy. 

The difference |V;;—V,';| is not so simply shown to be small. Indeed, this 
difference yields the summation which appears in the growth condition G,. To 
begin with, — Cr p| =| (Fee, Yes) —E(Ev pg, Vera)|, SO that, since » >No, 


| Cup Ty Cy p| = [7 (E:,> Vig) — HEws> Vira) | + 2Ec 


| (E.,66, Vig, Veo Vig) — (Ex, &5, Vep, Ps Vea)| 

+ |# (Ee &B, View Veo, Ver) — (Se a Via» Ve» Ver)| 

=“ | f(& i 1&8) Vip, Vip, Vie) — ft (ER »&, Ve 'p» Vip Vee) | 

+ |f(Ee,E3, Vir, Virw, Ves) — (Ee. €8, Vir, Virw, Virw)| + 2 


or 


| Cag — Crp] S41 + Ant As + Agt 2€c; 


where in the last inequality the A; are the respective terms in the previous 
sum. A bound will now be found for each 4;. 
From the comment following (4.2), and the uniform continuity of /, each 
of A,, A, and A, can be made less than 4¢, by taking | &; —&}|<d,. 
By the Lipschitz property of /, 
A, <L|Vig— Viral ) 
and thus one finds that 


(4.3) | Cee — Cve| Sat 2ect+L|Vig— Viral. 
Hence, using also the continuity of o’, one finds from (4.2) that, given ¢> 0, 


aol 
| Ve Vill (fi 2ec) a+ E 2 \Vee = Vera| Ai. 


provided that |&—&,| is less than an appropriate 6, and »>Nc. Because of 
the uniform continuity of f, o, 0’, t and 1’, this holds throughout ”. 


6* 
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Letting €,,=e+ (&+2€c) a*, one has 
es ; 
|| Ves — Verw|| SS enw + P| Vip — Vera ll Ap 
so that V satisfies the growth condition G,. 
Similarly, for suitable 6, and No, if | é? —é7| <0, and y > No, then ||Vi,, — Viv |S 


kl 
ew +L » |\Var—Var 
a=0 
Further, the hypothesis (d) of Theorem 2 does hold because of the uniform 
continuity of /, a, o’, t and t’. This completes the proof of the lemma. 


A} so that V satisfies the growth condition Gy. 


Theorem 2 can now be applied to the sequence of vines {V”} supplied by 
Lemma 2. Thus, there is a subsequence {V”} of {V”} which converges uniformly 
as 7->0o to a continuous function V: ¥~¥W. Because Y is (here) three-dimen- 
sional real space, the function V defines three real-valued functions V°, V1 and 
V2, the components of V, which are again continuous and are the uniform limits 
of V°”, V'” and V2”, respectively. The next lemma asserts V® is a solution of 
2-HP. 

For this remaining lemma, the more common (x, y)-notation for two inde- 
pendent variables will be used since it is somewhat more convenient for the 
proof, and because this lemma provides the final connection of the preceding 
theory with 2-HP. Thus, % will be given by {(x, y): OS “Sa, OS yS b\, while 
M’ will be given by the lines x=&%, y=nj. Also, f(x, y, V(x, y)) will mean 
(x,y, V(x, 9), V(x, »), V2(x, 9) 

Lemma 3. The components V°, V1, V? of V form a solution of the integral 
equations (1.7). 


Proof. Let 


ze 3; 


u(x, ¥) =o(x) +t(y)—o(0) +f St(s,é V(s,t))dsdt, 


0 0 


p(x, y) =o"(x) + ie (x, t, V(x, t)) dt, 
q(x, y) =t'(y) +f (s, y, V(s, y)) ds. 


Note that «, p and g are continuous (and therefore uniformly continuous) on &. 
To prove the lemma it must be shown that V°=u, Vi=p and V?=g onall of &. 
Using the equation defining V®”, that is, the first of equations (4.1), one finds 


gh ny 
[Ver — 425591) =e p1&, An; — jaPacenaee p) dsdi| 
0 0 
koi tt ae "B41 Std "B44 
“ a ei J Copdsdt— f ff f(s, t, V(s,t))dsdt 
“=o Folie af a 


—1 /-1 Sel MB 44 


Sy C5 -Sh(S ay 
S24 oe p—f(s,t,V(s,2))|ds de. 


VA 
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Because / and V are uniformly continuous on %, given e>0, the mesh can 
be made fine enough, say y>N,, so that (for Ss< Gar, 4p St<g41) ify >N,, 
then |Cig—f(s, #, V(s, t))|S|Cye—H(E., np, V,,)|+-¢. But the V” converge uni- 
formly to V. Thus, if »>.N,, then, again because of the continuity of f, 

| Cup ae / (s, t, V(s, t))| SS | Cup Ant ice NB» Vis) oe 2€ 


which can be written, using the definition of Cup» as 


|@"(:.78» Vas) — £ (E18 Vap)| + 2¢. 


Further, because @ converges uniformly to /, there is an N3 such that if »> Ng, 
then 


(4.4) | (&. np» Vas) — (Ei. 78, Ves) | <e. 
Hence 
pike jal 
[Vat — u(Sm)| Sd D3ed& Ang S3ead. 
a=0 B=0 


Since V°” converges uniformly to V° it follows that if y>N,, then 


| Vix — 4 (Es, m1)| S|Var— Ver'| + [Var — 4 (Ei n)| Set+3ead. 


Finally, because V° and w are uniformly continuous, if the mesh.@’ is fine enough, 
say vy >N,, so that for every (x,y) in 2 there is a node (&, 77) for which 
l| (x, ¥) —(&, 77) || is small, then 


|V°(x, y) — u(x, v)| S3e+3ead, 


and hence V°=u on all of 2. 


Precisely parallel arguments hold for V! and V?, so that Vi=# and V?=q 
on all of 2. This completes the proof of the lemma. 


For applications in the next section and for the Runge-Kutta procedure 
mentioned in the Introduction, the following is important. 


Remark. In Lemma 2, for each y, only a single function @’ (x, y, v) was used 
to define all three of V°”, V*” and V®” (see eqn. (4.1)). In the proofs of Lemmas 2 
and 3 the fundamental property of the sequence {@’} was that of its uniform 
convergence to f (cf. the reasoning employed between (4.2) and (4.3) in Lemma 2, 
and the inequality (4.4) in Lemma3). Instead of the single sequence, three 
sequences {6°"', {@'"! and {@°”} could just as well be used, so that in equations 
(4.1) C” would appear in the summation for V’”. That the three functions V” 
are well defined is easily checked, as before. The important requirement on the 
sequences {@’} (other than that each function is, of course, constant on each 
cell of @*”) is that they converge uniformly to fon 2 x¥. The key equations 
(4.3) and (4.4) then still hold. 

The result of the preceding Lemmas 2 and 3, and the remark just made, 
may be summed up in the following theorem. 


Theorem 3. Let % be the set {(x, y): OS “Sa, OS ySX0b} im the real plane; 
let V be 3-dimensional real space with points (v°, v!, v2). Let f be a real-valued 
function defined on & XW such that f is bounded and continuous and satisfies 
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the Lipschitz condition, with Lipschitz constant L, 
Lieu, Woy o) = [Ve v)| < L {|v — 7| +|v?— v?|} 


on all of LXV. Let o and v be continuously differentiable real-valued functions 
defined on 0<x<a and OS yXb respectively, and suppose o(0)=T (0). 

Let {M" and {.M*"} be convergent sequences of meshes in & and 4 XV, re- 
spectively, and {@\ (;=0, 1, 2) sequences of veal-valued functions defined on X XW 
such that €” is constant on each cell of M*”, and such that {6"} converges uniformly 
to fas yoo (7 fixed). 

Let N” denote the set of nodes of @’. 

Then, a sequence of functions V’: V”-—-¥V 1s defined by equations (4.1) with 
the property that a subsequence converges uniformly to a continuous function 
V: £+V. The first component V° of V is a solution of the differential equation 


Uny =f (%, YU, Uz, Uy) 
on & and satisfies the characteristic initial conditions 
u(x,0)=a(x) and u(0,y)=T(y). 


The second and third components V1 and V2 are, respectively, the x and y derivatives 
of the solution V°. 
§ 5. Convergence to a Given Solution 


It is well known that if the Lipschitz requirement on / in Theorem 3 is re- 
placed by the stronger requirement that 


lee Rue ee, Ojo (x, 00h my) Al a? | +|vt— o4| + |v24 =i 


then the solution to 2-HP is unique, while under the weaker condition in Theorem 3 
there may be several solutions. (A simple illustration is given, for example, 
by D14z [3], p. 361.) It is of interest to note that every continuously differentiable 
solution can be obtained as the limit of a sequence of approximate solutions of 
the type described in the preceding sections. This was shown by ZwIRNER [10] 
for the case when f=/(x, y, u). i 


Theorem 4. Let 2, Y, f, o and t be as in the hypothesis of Theorem}. Let 
u(x, ¥) be an arbitrary solution on & of 


(5.1) Uny =1(%,Y,U, Uz, Uy), U(x,0)=o(x), w(0,y)=T(y), 
where u, u, and u, are continuous on &. 


Then there is a sequence of vines of the type considered in Lemma 2 such that 
this sequence converges uniformly in & to a function V: ®>+¥W, and the first 
component of V is u(x, y). 


Proof. Let {4} and {.@*"| be convergent sequences of meshes on # and 
2% x’ respectively, defined by x=&, y=nj, vi=a!? (7=0, 1, 2). In general 
terms, all that is needed to prove the theorem is, for each », to define “suitable” 
functions @°”, @!” and @?” to be used respectively in the three equations (4.1). 


To this end, let p=u,, g=ty, Mapu (E.,15), Pap =P (EL NB)» dup = (Em) 
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and define 


(5.2) Ong Sg ; 
os Cl” = Put1,8 — Pap 

(5.3) ipa Pict Pet, 

5.4 C2, = dab Vap 

(5.4) — 


To define @”, on each cell of -@*” containing points (x, y,u (x,y), P(x, ¥),¢(x,¥)) 
for which €,S*< 1 and 43S y< 541, let the value of @” be C%,. On each 
remaining cell of @*” let the value of @” be any value taken on by io that cell. 

From the ordinary mean value theorem of ee calculus, it is evident that 


v 0 Ve say 
Cub = Gy P o> 1B) and Cee = q( em) 


np <p << NB+1 and a4 << gE < byes 
From a similar mean value theorem for functions of two independent variables 
(see Goursat [4]) it follows that 


Ov oe ) Z ‘ i 
Crtaetoy ayes Sines where BE ae os Ne <p” <NB41- 


Thus, using the differential equation (5.1) (note that £,=¢,—w,,), one finds 


Coy =F (EE nf, CER nf), PCER mB). (ER? PP), 
Cun = 1 (Ei Tp (Ean), (Eas TB). (Ea, Tp), 
Cop 1 (Ei mis, «(5.m5), B(E. 8)» (E18) 
From the continuity of f, uw, p and q it is thus clear that the sequences {@”} 
(7=0, 1, 2) converge to jas » oo. 

It remains only to be shown that the functions V°”, V'” and V?” defined 
by equations (4.1) using the above functions @°’, @'” and @®” respectively, do 
converge to u, p and q. In fact, it will be shown that V9” =u (&%, 11), Ver’ =P (Ek, 11) 
and V2”=4q(&&, 7). The last two of these equalities are established immediately 
upon the substitution of the formulas (5.3) and (5.4) for Cz and CZ into the 
last two of equations (4.1) respectively. It is a little more difficult to derive 
the first equality. Formula (5.2) for C}% substituted in the first of equations 
(4.1) yields the result (in the present notation) 


where 


eee 
Ov \) 
(5:5) Var = 0, +%—Oo+ 2, Dd (Ha41,641 — Ya41,6 — Ma,p+1 + Map) - 
o=—0 6—0 
In this, 
pai) Peal oer! jaa, (eal (Na al 
pa Uy +1,B+1 — Dimi Un p= » Uy Bi Di Mea Ung r Unis 
a=0 B=0 a=1 p=1 p=1 cL pe 
aoa eye pi k—-1 pat 
7 Se 
> ic 13 eae Dns Uy p= 2 Ugg 2 4e0 7 Li Ung Yeo» 
=O p—0 C=, p=O g=l gel Ch p= 
jaa pat alan ee Boil faa k—-1 Tea 
— Th 
2 Che a = eps Uy p= peed Ue Bi Ane Up Uo1> 
a—0 B=0 a=0 p=1 cal (a al ae 
k-1 1-1 k—1 1-1 1— 
Uy B aT Me Uy B tant 2 Hop + Uo0- 
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expressions on the right may now be substituted for the appropriate sum- 


mations in (5.5). All the terms but #,, then cancel on the right of (5.5) since 
Ung = Op, Upp = T, aNd My) =O. Thus Vy)’ = u,; as asserted. 
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